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© Method of designing three dimensional electrical circuits. 



© A computer aided design package is used to 
create a mathematical representation of a three- 
dimensional object. This object is defined as a set of 
surfaces oriented in space. A map of the flattened 
object is created by concatenating selected ones of 
the surfaces on a single plane. The outline of this 
map is then used in a computer aided circuit layout 
package as a printed circuit board on which an 
electrical circuit is placed and routed. The circuit is 
translated into a three-dimensional form correspond- 
ing to the surface of the object by translating and 
rotating the representation of the object to align each 



selected surface with the circuit description gen- 
erated by the circuit layout package. The portion of 
the circuit corresponding to the surface is then trans- 
ferred to a three-dimensional data structure having a 
format that is compatible with numerically controlled 
machining apparatus. This data structure is used to 
drive a numerically controlled phototool which cre- 
ates a three-dimensional mask that may be used to 
print the circuit on the surface of the three-dimen- 
sional object. 
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METHOD OF DESIGNING THREE-DIMENSIONAL ELECTRICAL CIRCUITS 



The present invention relates to the design of 
electrical circuits, and specifically to a method of 
producing electrical circuits on the surface of three- 
dimensional objects. 

Many electrical or electronic circuits are made 
by means of printed circuit boards (PCBs) In which 
electrical conduction paths are defined as traces of 
copper on the surface of a flat sheet of insulating 
material. In mire complex versions, there may be 
multiple layers of circuit, for example,, on the two 
sides of the PCB, which are connected to each 
other by plated through holes or vias in the PCB. 

PCBs provide both electrical connection and 
mechanical support for electronic components such 
as resistors, capacitors and integrated circuits, in 
addition, PCBs may include connectors used to 
couple the circuit, for example, to power sources, 
transducer or other circuits. Generally, PCBs are 
protected and supported in cases or other struc- 
tures to make usable appliances. 

In contradistinction to the flat, two-sided nature 
of the conventional PCB, a molded circuit board 
(MCB) is generally a three-dimensional object, hav- 
ing multiple surfaces. Thus, an MCB may combine 
both mechanical features and electrical functions 
into one artifact. 

This combined functionality is not achieved 
without some cost, since the design of such an 
MCB combines the mechanical design technology 
of molded or solid parts with circuit layout methods 
used to make PCBs. While there are readily avail- 
able computer software packages of program 
suites than can aid designers in either of these 
fields, these programs are generally incompatible 
and, thus, difficult to use in combination. This in- 
compatibility is fundamental to the design meth- 
odologies and is not merely an issue of translating 
between different data formats. 

The software packages for circuit design allow 
information about the placement of components in 
the PCB and the required connections among 
these components to be provided as an input data 
file, such as a netlist. The subsequent work of 
laying out the traces and vias is contained by the 
obligation to meet topological restrictions imposed 
by the input data file. Many of these program 
packages provide automatic trace routing algo- 
rithms which do some or all of the job with varying 
degrees of competence. All packages allow manual 
routing, in which traces and vias may be moved on 
the PCB while the circuit connections are pre- 
served. Neither the automatic algorithms nor the 
ability to move traces ad vias is a part of the 
available methodology of conventional mechanical 
engineering design packages used to design three- 



dimensional objects. Conversely, the circuit layout 
and design packages do not admit of the existence 
of a third dimension, except in that they may 
provide for several parallel layers of circuit traces 

5 and the local connection between them. 

The interrelationship between electrical and 
mechanical design s is very strong. Even in the 
design of flat circuit boards, the through holes 
which link circuit traces on opposite sides of the 

70 board and the holes which are used to support or 
attache components to the board are recognized as 
existing on all layers. Furthermore, many of these 
holes (i.e., vias) are logically associated with the 
traces. 

75 When a designer adjusts a design in the elec- 

trical design package by moving the trace and, 
thus, causing a via to be moved, the computer 
software attempts to preserve the electrical-me- 
chanical relationship by adjusting the routing on all 

20 layers of the board to conform to the new position 
of the via. Thus, in the electrical packages, a me- 
chanical feature is automatically adapted to circuit 
requirements. In the mechanical engineering de- 
sign packages, such topological relationships are 

25 not recognized. Mechanical features such as pat- 
terns of holes may be moved as a group, but in 
general, there is no way of relating these features 
to the continuity of adaptable decorative details on 
the surface of the object 

30 In attempting to extend circuit layout to sur- 

faces in three dimensions, a new set of problems 
arise due to the limitations of the two-dimensional 
circuit layout tools. When a solid part is to have 
circuit traces drawn, for example, on its inner and 

35 outer surfaces, electrically identical traces on the 
two surfaces may be very different due to the 
thickness of the part near edges or corners, or due 
to features that break the surface on one side but 
are not replicated on the other side. Thus, features 

40 such as stiffening ribs, standoffs, outer walls and 
mounting structures may make the topologies of 
two sides of an object very different. 

In the implementation of the design of a mol- 
ded part, tools are generally used which carry, in 

45 mechanical form, the image of what would be the 
artwork (i.e., traces and connecting pads) in normal 
two-dimensional circuit processes. In some cases, 
these tools may be generated using normal circuit 
board processes, for example, when the part in- 

50 eludes flat areas which can be dealt with as simple 
locally flat boards. However, to obtain full benefit of 
integrating mechanical and electrical designs, it is 
desirable to allow the traces to follow shaped sur- 
faces. To do this, the tools which are used to form 
the part should impose the traces on the part. 
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However, these tools are themselves solid ob- 
jects which ate fabricated by mechanical means 
such as machining or casting. Consequently, the 
traces and vias which constitute the electrical cir- 
cuit should be specified in whatever system is 
used to fabricate the tools, be it mechanical draw- 
ings interpreted by a skilled machinist or input data 
files which are used to direct the motion of a 
numerically controlled machine tool. 

The present invention is directed to a method 
for generating, from a representation of a two- 
dimensional circuit, a representation of a three- 
dimensional circuit which conforms to selected sur- 
faces of a representation of a three-dimensional 
object, comprising the steps of 

a) aligning one surface of the representation of 
the three-dimensional object in a coplanar rela- 
tionship with the representation of the two-di- 
mensional circuit; 

b) establishing a correspondence between said 
one surface of the three-dimensional object and 
a portion of the representation of said two-di- 
mensional circuit; 

c) associating circuit features of said portion of 
the representation of the two-dimensional circuit 
with said one surface; and 

d) copying said associated circuit features to the 
three-dimensional representation of said circuit. 
Figure 1 is a perspective view of an exemplary 
three-dimensional object suitable for use with 
the present invention. 

Figure 2 is plan drawing which illustrates the 
unfolding of the object of Fig. 1 to create a two- 
dimensional image of the object. 
Figure 3 is a circuit trace diagram placed and 
routed onto the image shown in Fig. 2. 
Figure 4 is a perspective view of the object 
shown in Fig. 1 with the circuit represented by 
the circuit trace diagram of Fig. 3 mapped onto 
its surface. 

Figures 5 to 8 are flow-chart diagrams which are 
useful for describing the present invention. 
The design process of an MCB begins with the 
creation of a three-dimensional model of the part. 
For the described embodiment of the invention, 
this three-dimensional model is generated in an 
Apollo workstation using the GMS CAD-CAM sys- 
tem of Unisys Corporation as the mechanical de- 
sign package. In this system, a three-dimensional 
object is specified as a set of quadrilateral surfaces 
in three dimensions. These surfaces are also 
known as ruled surfaces. The mechanical package, 
produces a data file describing the three-dimen- 
sional object. In the present embodiment of the 
invention, this file is in Initial Graphics Exchange 
Specification (IGES) format. An exemplary object is 
described below with reference to Fig. 1. 

The next step in the process is to "unfold" this 



three-dimensional object to generate a two-dimen- 
sional map of the ruled surface which, in the final 
artifact, are to include circuitry. In the present in- 
vention, this step is performed by concatenating 

5 the ruled surface, generated using the GMS pack- 
age, to produce an outline representing a two- 
dimensional PCB. As set forth below, although this 
outline should approximate the outline of the con- 
catenated ruled surfaces, it is not necessarily 

w isomorphic with the three-dimensional object. The 
process of generating an outline of this type is 
described bellow with reference to Fig. 2. 

The two-dimensional outline is generated, from 
data describing the concatenated surfaces, in an 

15 electrical design package, such as the Visula sys- 
tem produced by Racal Redac. This design pack- 
age is also used to specify the placement of fea- 
tures such as connecting pads and vias which 
define both electrical and mechanical coupling of 

20 circuit components. When these features have 
been placed, the electrical package is used to 
define a routing for traces that connect these pads 
and vias as indicated by the circuit description. For 
this embodiment of the invention, the placement is 

25 accomplished by an operator specifying the posi- 
tion of vias and the position and size of pads on 
the two-dimensional surface. The pads are iden- 
tified as being associated with specific components 
or with specific networks in the netlist. For the 

30 exemplary electrical design package set forth 
above, the routing operation is performed automati- 
cally to link the connecting pads and vias as in- 
dicated by the netlist file. The electrical package 
produces a data file describing the placed and 

35 routed circuitry. In this embodiment of the inven- 
tion, this file is a GERBER photoplotter file. Fig. 3 
is an example of a circuit routing produced by the 
electrical package. 

The GERBER file describes a connecting pad 

40 as a point on the two-dimensional surface and an 
aperture for a beam of light which is flashed on to 
produce an image of the pad on a photosensitive 
material. A trace is defined as a starting point, an 
ending point and an aperture for a beam of light 

45 which follows a path between the starting and end- 
ing points to produce an image of the trace. 

In the next step, the data files representing the 
three-dimensional object and the two-dimensional 
circuit diagram are combined by a program, which 

50 also runs on the Apollo workstation, to generate a 
data file representing a three-dimensional object 
having a three-dimensional circuit imposed on its 
surface. This program combines the two files by 
using space transformation algorithms which 

55 "rotate" a mathematical representation of the ob- 
ject in space to align each surface with a portion of 
the circuit diagram. As each surface is aligned, the 
circuitry circumscribed by that surface is removed 
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from the data file for the circuit diagram and added 
to a data file for the three-dimensional circuit. In 
the exemplary embodiment of the invention, this 
program produces an IGES file which describes a 
mask that may be applied to the object generated 
in the first step to impose the circuitry on that 
object, Fig. 4 is an illustration of an exemplary 
mask. 

In Fig. 1, the three-dimensional object shown in 
a perspective view includes seven surfaces: 101, a 
hidden surface as shown In Fig. 1, that is bounded 
by vertices A, C, E, N and O; 102, bounded by 
vertices A, B, J, G and C; 103, bounded by verti- 
ces C, E, F and H; 104, bouTided by vertices E, N. 
M and F; 105, bounded by vertices F, G, J, Q, P, 
M and L; 106, bounded by vertices B, R, Q and J; 
and 107, IThidden surface as shown in Fig. 1, 
bounded by vertices H, G, F and S. 

As set forth above, this representation of the 
object may be generated using the mechanical 
design package. The object is generated one sur- 
face at a time, where each surface is a quadrilat- 
eral. To specify a surface, two line segments are 
described in terms of X, Y and Z coordinates and 
the endpoints of these segments are joined to form 
the quadrilateral, it is noted that some of the sur- 
faces shown in Fig, 1 are actually composed of two 
or more quadrilaterals. For example, surfaces 102 
and 1 05 are composed of two and three 
quadrilaterals, respectively. 

The next step in the process is to generate a 
two-dimensional circuit layout area by concatenat- 
ing the surfaces of the object shown in Fig. 1. This 
process resembles the generation of a two-dimen- 
sional map pf the three-dimensional object. This 
step is performed using the mechanical package. 
Each surface of the three-dimensional object is 
copied and references to a common plane. The 
various surfaces are then oriented so that adjacent 
surfaces are joined. However, some surfaces which 
are contiguous in the three-dimensional object, are 
separated as a result of the mapping process. An 
exemplary map is shown in Fig. 2. 

A set of dimensions which describe the outer 
boundaries of the concatenated surface is applied 
as input values to the electrical package to create 
an equivalent two-dimensional map upon which the 
circuit is to be placed and routed. In order to link 
the map to the three-dimensional object one vertex 
is selected as a common orienting point. In this 
exemplary map, the orienting point is the vertex O. 

The map used by the electrical package does 
not need to be an accurate representation of every 
surface of the object. In Fig. 2, for example, the 
quadrilateral defined by the vertices D, J, G and H 
has been omitted from the map, and the surface 
107, defined by the vertices H, 5, F and G, has 
been incorporated in the map by lengthening the 



horizontal dimensions of the surfaces 103, 104, 
105a t 105b and 105c. 

The surface 107 is handled in this manner so 
that the map may have a regular structure suitable 

s for use by the electronics package. The inventor 
has determined that the routing achieved with this 
outline is more effective than that which may be 
achieved from a more isomorphically mapped out- 
line. Since the width of surface 107 has been 

70 accounted for by lengthening the surfaces 103 and 
105a, any circuit laid out in the portion of the map 
near the junction of these two surfaces will be 
transferred onto the surface 107. 

Circuitry may or may not be transferred to the 

75 omitted surface defined by the vertices D, J, G and 
H. This depends on the order in which the surfaces 
are evaluated when the two-dimensional circuit is 
applied to the three-dimensional object. However, 
the portions of an object which are to be occupied 

20 by circuitry may be limited to selected adjacent 
surfaces by concatenating only those surfaces to 
generate the outline. A vertex of one of these 
surfaces should then be chosen as the orienting 
point of the outline. 

25 For the purposes of the electrical package, the 

area onto which the circuit may be placed is repre- 
sented by the outline of the map shown in Fig. 2, 
with the respective pairs of vertices C and C and 
H and G' joined as indicated. The boundaries of 

30 the individual surfaces which were concatenated to 
generate this map are ignored by the electrical 
package. 

The circuit produces by the electrical package 
for this example is shown in Fig. 3. This circuit was 

35 generated by defining connecting pads and vias 
corresponding to two integrated circuits 302 and 
304 and a transistor 306 as well as various other 
pads and vias as indicated. In the electrical pack- 
age, the individual vias are labeled to correspond 

40 to networks, or nets, of a circuit defined by a netlist 
data file. The electrical package then routes the 
circuit defined by this netlist onto the outline gen- 
erated by the mapping process. Fig. 3 shows the 
result of this routing operation. 

45 In the circuit diagram shown in Fig. 3, traces 

have been placed in the void region defined by the 
vertices N, E, C, c\ e' and hT. If this were a flat 
board, traces such as these would generate an 
error indication from the electrical package. In this 

so instance, however, the edge defined by the vertices 
n', e' and c' is joined to the edge defined by the 
vertices N. E and C in the three-dimensional object 
and, thus, the void does not actually exist. Depend- 
ing on the electrical package used, some error 

55 checking may need to be deimplemented to allow 
traces to be routed across voids of this type. 

Care should be exercised where traces have 
been placed across voids, to ensure that the ends 
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of traces on opposite sides of the void are properly 
aligned for the three-dimensionai object For exam- 
ple, consider trace 308 generated by the electrical 
package. If the portions of this trace which extend 
across the void were omitted from the three-dimen- 
sional object, there would be no electrical connec- 
tion because the end of the trace at the edge N -E 
would not meet the end of the corresponding trace 
at the edge N-E. To correct this problem, the 
operator of the electrical package moves trace 308 
to the position 310 indicated in phantom. Other 
traces which extend across voids should be exam- 
ined for similar adjustment. 

If circuitry is to be placed on both sides of the 
three-dimensional object, each of the steps set 
forth above would be performed separately for the 
inner and outer surfaces of the object. In the 
present embodiment of the invention, the operator 
is responsible for ensuring that any differences in 
surface area due to curvature of the three-dimen- 
sional object or due to the thickness of its walls are 
properly represented by the generated circuit out- 
lines. 

Figures 5A and 5B are a flow-chart diagram 
which illustrates the overall process of creating a 
MCB for a three-dimensional part. Steps 502, 504 
and 506 relate to the design of the three-dimen- 
sional part, using the mechanical design package, 
and the creation of the IGES file describing the 
part. Steps 508 and 510 concern the creation of the 
two-dimensional surface outline from the three-di- 
mensional part description. Steps 512, 514 and 516 
relate to processes performed using the electrical 
package. These steps in Figures 5A and 5B are 
described in detail above. 

Step 518 of Fig. 5B is implemented, in the 
present embodiment of the invention, by a program 
written in the computer language SMALLTALK. A 
source listing of this program is included as an 
appendix to the present application. The step 518 
is briefly described above and is described in 
detail below with reference to Figures 6A to 8. As 
set forth above, this step produces an output data 
file which describes the circuit in three-dimensional 
space as it has been mapped onto the three- 
dimensional part. 

In step 520, a test is made to determine if the 
output file is to be used to drive imaging equip- 
ment, in which case the file is sent to the imager at 
step 522. Otherwise, at step 524, the file is used by 
the mechanical design system, together with the 
IGES file that describes the three-dimensional part 
(provided as an input 526 to the step 524), to 
generate, at step 528, commands for a phototool. 
The phototool, in turn, generates a three-dimen- 
sional mask which may be used, as set forth 
above, to print the circuit on the surface of the 
three-dimensional part. 



Figures 6A and 6B are a flow-chart diagram 
which provides grater detail on the process of 
mapping the two-dimensional circuit onto the three- 
dimensional object (step 518). 

5 The first step in this program, 602, is to create 

an aperture list. This list is generated by the oper- 
ator from data used to produce the GERBER circuit 
description file in the electrical package. This list 
defines the sizes of light beam apertures to be 

w used by the photo-plotter to create images of the 
pads and traces on the photosensitive output me- 
dium. 

In the next step, 604, the operator creates a 
surface list, defining the sequence in which sur- 

15 faces of the three-dimensional object are to be 
evaluated to receive circuitry from the two-dimen- 
sional circuit diagram. This list is created by the 
operator from a list of surface identifiers for the 
three-dimensional object. These identifiers are tak- 

20 en from the IGES file. 

At step 606, the program reads in the data 
describing the circuit from the GERBER data file 
606. This data is put into three-dimensional format 
at step 610. In this step, each pad and each trace 

25 are defined as respective ruled surfaces, if a format 
compatible with the IGES file but specifying only 
two dimensions. The third dimension is provided 
later in the procedure. This step also generates 
separate TRACE and FLASH lists which include 

30 arrays of records which define the positions and 
dimensions of the respective traces and connecting 
pads in the two-dimensional circuit. A separate pair 
of lists is established for each surface of the three- 
dimensional object, however, at this time, all lists 

35 are empty except for the respective FLASH and 
TRACE lists that correspond to the first surface in 
the surface list. These two lists are initialized to 
include all of the respective flash and trace records 
for the two-dimensional circuit. 

40 The next step, 612, reads the mechanical sur- 

face data from the IGES part file 614. In step 616, 
reached via the offpage connector 2, the three- 
dimensional circuit representation is mapped onto 
the three-dimensional part using the data from the 

45 IGES file. This step adds the third dimension to the 
circuit description that was generated at step 610. 

Step 618 illustrates the possible output data 
formats of the three-dimensional circuit description 
generated at step 616. In the present embodiment 

so of the invention, the output data is in the format of 
an IGES file. However, as illustrated in Fig. 6B, 
output data files may be prepared in several two- 
dimensional and three-dimensional formats. 

Figures 7 and 8 are flow-chart diagrams which 

55 show details of the steps 616, described above. 
Fig. 7 is a flow-chart diagram of a procedure MAP, 
which is invoked at step 616, and Fig. 8 is a flow- 
chart diagram of a procedure CLIP which is in- 
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voked from the procedure MAP as set forth below. 

The procedures MAP and CUP fit to two-di- 
mensional circuit onto the surface of the three- 
dimensional object. The fitting process may be 
analogized to the wrapping of the three-dimen- 
sional object with a sheet of paper representing the 
circuit. As each surface of the object is aligned with 
the sheet of paper, all circuitry which covers the 
surface is clipped out of the paper and saved The 
clipped circuitry is marked as to the surface to 
which it belongs. This marking process establishes 
a dependency between portions of the circuit and 
respective surfaces of the three-dimensional object, 
That is to say the location in space of a section of 
the circuitry is dependent on the position in space 
of the surface onto which it has been mapped. 
Following this analogy, the folding and clipping 
continues until there is not more circuitry to be 
mapped or until there are no more surfaces to be 
mapped onto. When this process is applied to a 
three-dimensional object having both interior and 
exterior surfaces, each of these surfaces is 
mapped separately and the results of the separate 
mappings are merged as a final step. 

The first step in the MAP procedure is to select 
a surface from the surface list. If this is the first 
invocation of the procedure, then first entry in sur- 
face list is selected. Otherwise, the seiected sur- 
face is determined by the invocation of the MAP 
procedure as set forth below. 

At step 704, the selected surface is examined 
to determine if it has been used, if so, or if the 
TRACE and FLASH lists for the surface are empty, 
step 706 is executed to select the next surface. 
The FLASH and TRACE lists for the surfaces were 
established at step 610 above. 

When a suitable surface has been selected, 
step 708 is executed to "rotate" and "translate" the 
selected surface into the X,Y plane. These oper- 
ations are a mathematical translation of the repre- 
sentations of the three-dimensional object to 
achieve a surface representation having a normal 
vector which is parallel with the Z-axis. An exem- 
plary method of performing the rotation is de- 
scribed in chapter 22 of a textbook by W.M. New- 
man et al, "Principles of Interactive Computer 
Graphics", McGraw Hill, 1979, pp. 333-354. 

During the rotation and translation operations, 
the orienting point of the three-dimensional object 
is maintained as a reference point for the two- 
dimensional circuit diagram. Thus, when these op- 
erations are complete, there is a one-to-one cor- 
respondence between the target surface and a 
portion of the circuitry. 

In step 710, the portion of the circuitry which is 
aligned with the target surface is clipped from the 
data structure created in step 6120 and made 
dependent on the target surface. The clipping op- 



eration is described in detail below, in reference to 
Fig. 8. 

In the next four steps, 712, 714, 716 and 718 
the MAP procedure invokes itself recursively for 

5 each surface which is adjacent to the target sur- 
face, as indicated by the surface list. The steps 
712 to 718 are executed if the surface list indicates 
that the target surface has respective left, right top 
or bottom adjacent surfaces. These recursive in- 

w vocations ensure that all surfaces of the three- 
dimensional object will be evaluated in the map- 
ping process. 

When all surfaces have been evaluated or 
when no remaining surface has entries in its 

is TRACE or FLASH lists, step 720 terminates the 
mapping and clipping processes. 

Figure 8 illustrates the clipping procedure in- 
voked in step 710 above. The clipping algorithm 
used in this embodiment of the invention is gen- 

20 eraily the same as the Cohen-Sutherland algorithm 
which is set forth in chapter 6, pp. 65-67 of the 
above-referenced text by W.H. Newman et al. 

In the first step, 802, the left, right, top and 
bottom surfaces, if they exist, which are adjacent to 

25 the selected surface are identified. Next, at step 
804, a trace r is selected from the TRACE list, (T- 
[s]) for the selected surface s. At step 806, the 
procedure determines if the trace t is entirely with- 
in the clipping polygon defined by the selected 

30 surface. If so, a transformation matrix dependency 
for the trace r is set, at step 808, to indicate that 
the trace is mapped onto the surface s. Next, at 
step 810, the trace r is deleted from the TRACE 
list T[s] and placed in the master trace list. 

35 If the trace r is not entirely within the clipping 

polygon, it may either be entirely outside of, or 
partly inside and partly outside of the polygon. This 
determination is made at step 812. If the trace is 
entirely outside the polygon, step 814 determines 

40 which edge of the polygon is closest to the trace 
and, at step 816, assigns the trace to the TRACE 
list for the adjacent surface which shares that edge 
with the clipping polygon. 

If the trace t is only partly outside of the 

45 clipping polygon, the trace is clipped at the edge of 
the polygon to form two traces, r and t and both 
traces are inserted into the TRACE list T[s] for the 
selected surface s. 

If, after steps 810 or 816, step 818 determines 

50 that the last trace has not been processed; or 
unconditionally after step 822, step 823 is executed 
to select the next trace to be used as the trace r. 
Step 823 branches to step 806 which processes 
the trace r as set forth above. 

55 If, at step 818, it is determined that the TRACE 

list for the selected surface is empty, step 824 is 
executed which processes the FLASH lists for the 
surface in substantially the same manner as de- 
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scribed above. The processing of the FLASH lists 
is the same as the processing of the TRACE lists 
except that there is no analogue to steps 820 and 
822 since, in this embodiment of the invention, a 
flash (connecting pad) may not be split across two 
surfaces. It is contemplated, however, that the sys- 
tem may be modified to accommodate such a split. 

The final step, 826, in the clipping procedure is 
to mark the selected surface s as used. This mark- 
ing is tested in step 704 as set forth above. 

It may desirable to specify that some ruled 
surfaces of a three-dimensional object are not to be 
imprinted with circuitry. This may be accomplished, 
for example, by marking these surfaces as used 
before the mapping and clipping processes are 
performed. 

A method of producing a printed circuit on the 
surface of a three-dimensional object has been 
described. While this invention has been described 
in terms of a preferred embodiment, it is con- 
templated that it may be practiced as outlined 
above within the spirit and scope of the appended 
claims. 



f) repeating steps a) through e) until no more 
circuit feature of the representation if the two- 
dimensional circuit can be aligned with the se- 
lected surfaces of the representation of the 
5 three-dimensional object. 

5. The method of claim 1 wherein step a) is pre- 
ceded by the steps of generating the representa- 
tion of the three-dimensional object as a set of 
surfaces, including said selected surfaces, where 
w each sur face in said set of surfaces is defined in 
terms of a shape and a position in space; mapping 
the representation of the three-dimensional object 
onto a two-dimensional map by concatenating the 
shapes of each of said selected surfaces; placing 
15 circuit features corresponding to components of a 
two-dimensional circuit on said two-dimensional 
map; and routing circuit features corresponding to 
connections of components of said two-dimensional 
circuit on said two-dimensional map to generate 
20 said representation of the two-dimensional circuit. 



Claims 



25 



1. A method for generating, from a representation 
of a two-dimensional circuit, a representation of a 
three-dimensional circuit which conforms to select- 30 
ed surfaces of a representation of a three-dimen- 
sional object, comprising the steps of 

a) aligning one surface of the representation of 
the three-dimensional object in a coplanar rela- 
tionship with the representation of the two-di- 35 
mensional circuit; 

b) establishing a correspondence between said 
one surface of the three-dimensional object and 
a portion of the representation of said two-di- 
mensional circuit; 40 

c) associating circuit features of said portion of 
the representation of the two-dimensional circuit 
with said one surface; and 

d) copying said associated circuit features to the 
three-dimensional representation of said circuit. 45 

2. The method of claim 1 further including the step 
of repeating steps a) through d) for each selected 
surface in the representation of the three-dimen- 
sional object. 

3. The method of claim 1 wherein step a) is pre- 50 
ceded by a step of defining an order in which the 
selected surfaces are processed by steps a) 
through d). 

4. The method of claim 1 further including the 
steps of 55 

e) deleting said associated circuit features from 
the representation of the two-dimensional circuit; 
and 
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Model subclass: #Mcb 

instanceVariableNames: 'traces flashes geom range xrange yrange zrange transform slist mt mp ' 
class VariableNames: f * 
poolDictionaries: '* 
category: 'Graphics-Circuits'! 
Mcb comment: 

'This class will model a 3 dimensional 
molded circuit board and its traces/! 



!Mcb methodsFon 'accessing'! 

flashAdd: aPad 

"Add a pad to the collection of flashes" 

flashes add: aPad! 

flashes 

"Return the collection of flashes" 

t flashes! 

flashes: aCollection 

"Save the collection of flashes" 

flashes «• aCollection! 

geom 

"Return the collection of geometries" 

t geom! 

geomAdd: aGeom 

"Add geometry to the collection of geometries" 

geom <- aGeom! 

mp 

"Return the collection of mapped flashes (pads)" 

t mp! 

mp: aCollection 

"Save the collection of mapped flashes (pads)" 

mp *■ aCollection! 

mt 

"Return the collection of mapped traces" 

t mt! 

mt: aCollection 

"Save the collection of mapped traces" 
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nit aCollection! 
range 

"return the 3D point which is the absolute extents of the model" 
t range! 

range: a3DPoint 

"save the 3D point which is the absolute extents of the model" 

range «- a3DPoint! 

slist 

"Return the surface list which orders mapping" 

t slist! 

slist: aLinkedList 

"Save the linked list of surfaces" 

slist aLinkedList! 

traceAdd: aTrace 

"Add a trace to the collection of traces" 

traces add: aTrace! 

traces 

"Return the collection of traces" 

t traces! 

traces: aCollection 

"save the collection of traces" 

traces «- aCollection! ! 

!Mcb methodsFon 'initialize*! 

initialize 

"initialize a collection of traces and geometry" 

transform OrderedCollection new: 1. 
slist*- Array new: 10! ! 

!Mcb methodsFon 'operations'! 

change 

"modify a point on a trace" 

I aTrace aPoint bPointl 

bPoint «- Point3D x:0.5 y:0.0 z:0.5. 
aTrace «- traces at: 1. 
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aPoint <- aTrace beginPoint. 
aPoint «- aPoint + bPoint- 
aTrace beginPoint: aPoint, 
traces at; 1 put: aTrace. 
self changed! 

changed: aString 

"pass along the Gerber command line being processed" 



self changed: aString! 

getGerber 

"read in and process a Gerber file to 
initialize flashes and traces" 

I gerber I 

gerber «- Gerber new. 
gerber readFile. 

Transcript cr, show: ' Gerber file read complete \cr,endEntry. 
self flashes: gerber flashes, 
self traces: gerber traces, 
self range: gerber extents! 

getlges 

"read in and process a Iges file to 
ini tializ e surfaces" 

I iges I 

iges «- Iges new. 
iges readFile. 

Transcript cr, show: ' Iges file read complete * ;cr t endEntry. 
self geomAdd: iges surfaces! 

getSlist 

"get list of surfaces to process" 

I cm link slink aString aFileName aFile I 
cnt «- L 

link *• LinkedList new. 
slink «- SurfList new. 
(Binary Choice 

message: *Do you wish to read in a surface list?') 

ifTrue:[ 

aFileName «- FilllnTheBlank 

request: 'What is the name of the surface file?* 

initialAnswen 'surface.list' . 
aFile FiieStream oldFileNamed: aFileName. 
siist «■ (aFile fileln). 
aFile close.] 

ifFalse:[ 
[ BinaryChoice 

message: 'Create a surface list?'] 
whileTrue:[ [BinaryChoice 

message: 'Add a surface?*] 
whileTrue: [ aString «- FilllnTheBlank 
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request: 'Enter surface name/ 
initialAnswen 'SI', 
slink name: aString. 
link add: (slink deepCopy). 
slink release. 
aS firing release.], 
self slist at: cnt put: (link deepCopy). 
[link size > 0] 

whileTrue:[link removeFirst]. 

cnt*- cnt + 1J. 
(Binary Choice 

message: 'Do you wish to save the surface list?') 
ifTrue:[ 
aFileName release. 
aFileName «- FilllnTheBlank 

request 'Enter surface list file name for saving.' 

initialAnswen 'surface.list\ 
aFile «- FileStream fileNamed: aFileName. 
aFile reset 
slist storeOn: aFile. 
aFile close.]. 
]. 

Transcript cr; show: ' Surface list creation complete. ';cr,endEntry.! 

inspect 

"inspect the traces" 

traces inspect! 

mapping 

"map circuits onto surfaces" 

I aSurface mappedCircuits surfNamel 
surfName «- (self slist at: 1) first 
1 to: (geom size) do:[:cnd 

((self geom at: cnt) name « (surfName name)) ifTme:[ 
aSurface «- self geom at: cnt].]. 

mappedCircuits *- aSurface mapToNextSurface: self traces with: self flashes and: self geom and: self slist. 

self mt: (mappedCircuits at: 1). 

self mp: (mappedCircuits at: 2). 

Transcript cr,show: 'Mapping Comp3ete';cr,endEntry! ! 

It ^ W| 

Mcb class 

instance VariableNames: "! 



!Mcb class methodsFon 'instance creation'! 
new 

"initialize a collection of traces" 
I newMcb I 

newMcb «■ super new. 
newMcb initialize. 
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tnewMcb! ! 
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View subclass: #McbView 

instance VariableNames: 'volume dorigin projection * 

class VariableNames: " 

pooIDictionaries: M 

category: 'Graphics-Circuits'! 
McbView comment: 
'This class is for viesing a 3 dimensional 
molded circuit board/! 



IMcbView methodsFon 'displaying*! 

displayPad 

"Display the model's Pads" 

I setl aPad tVal aForm width dPad nPad hwl 

( self model mp = nil) ifFalse:[ setl self model mpj 

ifTrue:[setl «- self model flashes.], 
self dorigin: self insetDispiayBox bottomLefL "set origin to lower left corner of window" 
tVal «- setl size. 
1 to: tVal do: [ :index I 

aPad «- setl at: index. 

nPad «- aPad convToView: self volume, "convert world to display coordinates" 

dPad <- self projectPad: (self projection) with: nPad. "do projection calculations here" 

dPad «■ self scale: dPad. "scale the object to fit the window" 

width «• self findPadWidth: aPad. "set shape correctly for display" 

hw *• (width® width) / (2@2). " offset pad center 1/2 width" 

dPad <- dPad - hw. 

(aPad padshape = 'round') ifTrue:[ "if pad is round set size" 

aForm <- Form dotOfSize: width. 

aForm offset: dPacL] 
ifFalse:[ 

aForm «- Form new extent: width@width offset: dPad. "if pad is square set width" 

aForm black.]. " turn it black" 

aForm displayOn: Display 
at: self dorigin 

clippingBox: self insetDispiayBox.]! 

displays urf ace 

"Display the model's surfaces" 

! setl aSurface tVal aForm width cSurface bSurface hwl 

setl *• self model geom. 

self dorigin: self insetDispiayBox bottomLeft "set origin to lower left comer of window" 

tVal «- setl size. 

1 to: tVal do: [ :index I 

aSurface <- setl at: index. 

bSurface aSurface convToView: self volume, "convert world to display coordinates" 
"bSurface inspect" 

cSurface «- self projectSurface: (self projection) with: bSurface. "do projection calculations here" 
"cSurface inspect" 

cSurface «• self scale: cSurface. "scale the object to fit the window" 

aForm «- Form new extent: 1@ 1 . "make a form which could vary by width" 
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aForm black. "turn it black" 

cSurface form: aForm. "use the black form for display" 

cSurface displayOn: Display 
at; self dorigin 

clippingBox: self insetDisplay Box.]! 

displayTrace 

"Display the moders traces" 

I setl aTrace tVal aForm width aLine nTrace hwi 

( self model rat = nil ) LfFalse:[ setl «• self model mt.] 

ifTrue:[setl «■ self model traces.], 
self dorigin: self insetDisplayBox bottomLeft. "set origin to lower left corner of window" 
tVal *■ setl size. 
1 to: tVal do: [ :index I 

aTrace «- setl at: index. 

nTrace <- aTrace convToView: self volume, "convert world to display coordinates" 

aLine «- self projectTrace: (self projection) with: nTrace. "do projection calculations here" 

aLine «- self scale: aLine. "scale the object to fit the window" 

width *■ self findTraceWidth: aTrace. "set width correctly for display" 

hw «- (width@ width) / (2@2). "offset line points for form width" 

aLine beginPoint: (aLine beginPoint)- hw. 

aLine endPoint: (aLine endPoint)-hw. 

aForm «• Form new extent: width@ width, "make a form vary by width" 

aForm black. "turn it black" 

aLine form: aForm. "use the black form for display" 

aLine displayOn: Display 
at self dorigin 

clippingBox: self insetDisplayBox.]! 

displayView 

"Display the model's pads and traces" 

self volume: self model range, 
self displays urface. 
self dispIayPacL 
self displayTrace! 

redraw 

"redisplay upon command from controller" 
self display! 
update 

"redisplay upon change in model" 

self display! 

update: aParameter 

"Display the value of the model." 

I box pos I 
self clearlnside. 

"Position the text at the left side of the display area." 
box «- self insetDisplayBox. "get the view's box" 
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pos «• box origin + (4 @ (box extent y / 3)). 

"Concatenate the components of the output string and display them/ 
('command: \ aParameter) asDisplayText display At: pos! ! 



IMcbView methodsFor: ' initialization'! 



startup 

"set view variables to default values" 



self volume: 12@12@12. 
self dorigin: 0@0@0. 
self projection: I! ! 



IMcbView methodsFor. 'transformation 1 ! 



findPadWidth: aPad 

"set width for form so object is displayed at correct size" 



I tx wt nw I 

tx «- self volume x. 

wt <- aPad padsize. 

nw *• ((wt / tx) * 1000) truncated. 

nw < 1 ifTrue:[ nw <- 1J. 

:nw! 



findTrace Width: aTrace 

"set width for form so object is displayed at correct size" 



I tx wt nw I 

tx «- self volume x. 

wt <- aTrace w. 

nw «• ((wt / tx) * 1000) truncated, 
nw < 1 ifTrue:[ nw «- 1.]. 
tnw! 



projectPad: anlnteger with: aPad 

"cabinet projection for now. In future will control how the drawing is projected 
by changing the coordinates to match the required view" 



I nPad mat I 



mat «- Matrix new: 4 by: 4. 
mat setToZero. 
mat atPoint: 1@1 put: 1.0. 
mat atPoint: 2@2 put: 1.0. 
mat atPoint: 4@4 put: 1.0. 
mat atPoint 1@3 put: 0.4477. 
mat atPoint: 2@3 put 0.8941. 
aPad rotate: mat 

nPad «- aPad coordinates as2DPoint 
tnPad! 



projectSurface: anlnteger with: aSurface 

"dummy for now. In future will control how the drawing is projected 
by changing the coordinates to match the required view" 
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I mat I 

mat «- Matrix new: 4 by: 4. 

mat setToZero. 

matatPoint: 1@1 put: 1.0. 

mat atPoint: 2@2 put: 1.0. 

mat atPoint: 4@4 put: 1.0. 

mat atPoint: 1@3 put: 0.4471. 

mat atPoint; 2@3 put: 0.894L 

aSurface rotate: mat. 

aSurface el: (aSurface el as2DLine). 

aSurface e2: (aSurface e2 as2DLine). 

aSurface e3: (aSurface e3 as2DLine). 

aSurface e4: (aSurface e4 as2DLine). 

t aSurface! 

projectTrace: anlnteger with: a3DIine 

"cabinet projection for now. In future will control how the drawing is projected 
by changing the coordinates to match the required view" 

Imatl 

mat ♦* Matrix new: 4 by: 4. 
mat setToZero. 
matatPoint: 1@1 put: 1.0. 
mat atPoint: 2@2 put: 1.0. 
mat atPoint: 4@4 put: 1.0. 
mat atPoint: 1@3 put: 0,4471. 
mat atPoint: 2@3 put: 0.8941, 
a3Dline rotate: mat. 
ta3Dline as2DLine! 

scale: anObject 

"scale an object. It must understand the scaleBy: message" 

r anObject scaleBy: (self displayTransformation scale)! ! 

IMcbView methodsFon 'accessing'! 

dorigin 

"return the x,y,and z drawing origin of the display space in inches as a 3D point." 
t dorigin! 
dorigin: aPoint 

"set the x,y,and z drawing origin of the display space in inches as a 3D point." 
dorigin «- aPoint! 
projection 

"return the projection selector of the display space to be used in selecting function" 

tprojection! 

projection: anlnteger 

"return the projection selector of the display space to be used in selecting function" 
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projection «- anlnteger! 
volume 

"return the x,y,and z volume of the display space in inches as a 3D point." 
r volume! 
volume: aPoint 

"set the x,y f and z volume of the display space in inches as a 3D point." 
volume *■ aPoint! ! 



McbView class 

instanceVariableNames: "! 



IMcbView class methodsFor. 'instance creation'! 
open 

"Open a view for a new mcb." 
"McbView open." 

I mcbwView top View mcbw I 
mcbw Mcb new. 

mcbwView «• McbView new 
model: mcbw; 

controller. McbController new; 
borderWidth: 1; 
insideColon Form white. 

mcbwView startup. 

top View «- Standards ystem View new 
label: 'Mcb'; 
minimumSize: 200@200; 
addSubView: mcbwView. 

top View controller open! ! 
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MouseMenuController subclass: #McbController 

instanceVariableNames: " 

class VariableNames: " 

pooiDictionaries: " 

category: 'Graphics-Circuits'! 
McbController comment: 

'This class is the controller for constructing a 3 dimensional 
molded circuit board.*! 



IMcbControllermethodsFon 'menu messages'! 

change 

"tell model to change" 

self model change! 

gerber 

"tell model to read in a Gerber file." 
self model getGerber! 

iges 

"tell model to read in an Iges file," 

self model getlges! 

inspect 

"tell model to inspect traces" 

self model inspect! 

inspectview 

"inspect the view" 

self view perform: # inspect! 

mapping 

"tell model to map circuits." 

self model mapping! 

redraw 

"tell view to redisplay" 

self view redraw! 

slist 

"tell model to read in an surface list file." 

self model getSlist! ! 
! McbController methodsFon 'control defaults'! 
isControlActive 
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"Inherits control from superclass MenuMouseController 
As the blue button menu is used for reframing views only, 
mcb is activated by pressing the other two buttons" 

t super isControlActive & sensor blueButtonPressed not! ! 

IMcbController methodsFor: 'initialize*! 

initialize 

"McbController initialization" 

super initialize. 

self yeUowButtonMenu: (PopUpMenu labels: 
'Redraw 

Change 
Inspect 
Inspect View 
Gerber In 
IGES In 
Surface List 

Mapping') m 

yellowButtonMessages:#(redraw change inspect inspectview gerber lges slist mapping). ! 
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Object subclass: #Gerber 

instance VariableNames: 'g d xl x2 yl y2 oldApt size m apShape traces flashes apertures xmin xmax ymin 
ymax zmin zmax * 

class VariableNames: " 

poolDictionaries: " 

category: 'Graphics-Circuits'! 
Geiber comment: 

'This class processes Gerber photoplotter 
files and reconstructs the circuitry for 
use by class Mcb.'! 



IGerbermethodsFon 'parsing'! 

getDtaString 

"extract the D command from the command line 
and save it in d" 

I pString s I 
s *■ aString size. 
1 to: s do: [:index I 

(SD as (aString at: index) ) 

ifTrue:[pString «- String with:(aString at: index) with: (aString at: index + I) 
with: (aString at: index + 2). 
self d: pString.].]! 

getG: aString 

"extract the G command from the command line 
and save it in g" 

I pString s I 
s aString size. 
1 to: s do: [:index I 

(SG = (aString at: index) ) 

ifTnie: [pString *- String with:(aString at: index) with: (aString at: index + 1) 
with: (aString at: index + 2). 
self g: pString.].]. 

"self inspect."! 

getM: aString 

"extract the M command from the command line 
and save itinm" 

1 pString s i 
s*- aString size. 
1 to: s do: [:index I 

(SM = (aString at: index) ) 

ifTrue: [pString «- String with:(aString at: index) with: (aString at: index + 1) 
with: (aString au index 4- 2). 
self m: pString.].]! 

getX: aString 

"extract the X value from the command line 
and save it in x2, as the current x value. 
This method assumes the Gerber numbers 



20 



EP 0 416 568 A2 



are in 2.3 format, with or without leading zeros." 

I aFloat aNumber s aChar i i 
aNumber *• 0. 
s*- aString size. 
1 to: s do: [-.index I 

(SX = (aString at: index) ) 
ifTrue:[i*- index +1. 

[(aChar «• aString at: i) isDigit] whileTrue: [ 

aNumber «- aNumber * 10 + (aChar aslnteger - 48). 
i<-i+lj. 
aRoat «- aNumber / 1000.0. 
self x2: aFloat. 

( aFloat < xmin) ifTrue:[ self xmin: aFloat.]. 
( aRoat > xmax) ifTrue:[ self xmax: aFloat.]. 
].].! 

getY: aString 

"extract the Y value from the command line 
and save it in y2, as the current y value. 
This method assumes the Gerber numbers 
are in 2.3 format, with or without leading zeros." 

! aFloat aNumber s aChar i I 
aNumber «■ 0. 

aString size. 
1 to: s do: [:index I 

(SY = (aString at: index) ) 
ifTrue:[i*- index +1. 

[(aChar «- aString at: i) isDigit] whileTrue: [ 

aNumber «- aNumber * 10 + (aChar aslnteger - 48). 
i «- i+IJ. 
aFloat «- aNumber / 1000.0. 
self y2: aFloat 

( aFloat < ymin) ifTrue:[ self ymin: aRoat]. 
( aRoat > ymax) ifTrue:[ self ymax: aRoat]. 
]J! 

parse: aCommand 

"Take a command line and parse it into its proper components." 

self getG: aCommand. 
self getX: aCommand. 
self getY: aCommand. 
self getD: aCommand. 
self getM: aCommand! ! 

IGerber methodsFon 'accessing'! 

aperture 

"return the string which is the D aperture" 

tapenure! 

aperture: aString 

"save the string which is the D aperture" 
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aperture «- aString! 
apShape 

"return the string which is the aperture shape" 

tapShape! 

apShape: aString 

"save the string which is the aperture shape" 

apShape <- aString! 

d 

"return the string which is the D command" 
td! 
d: aString 

"save the suing which is the D command" 

d«- aString! 

flashAdd: aPad 

"Add a pad to the collection of flashes" 

flashes add: aPad! 

flashes 

"Return the collecuon of flashes" 

t flashes! 

flashes: aCollection 

"Save the collection of flashes" 

flashes «• aCollecdon! 

g 

"return the string which is the G command" 
tg! 
g: aString 

"save the string which is the G command" 
g «- aString! 

m 

"return the string which is the M command" 
tm! 
m: aString 

"save the string which is the M command" 
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m <-aString! 

size 

"return the number which is the aperture size" 

rsize! 

size: aFloat 

"save the number which is the aperture size" 

size «- aFloat! 

traceAdd: aTrace 

"Add a trace to the collection of traces" 

traces add: aTrace! 

traces 

"Return the collection of traces" 

1 traces! 

traces: aCollection 

"save the collection of traces" 

traces «- aCollection! 

xl 

"return the number which is the previous x location" 
txl! 
xl: aFloat 

"save the number which is the previous x location" 
xl «- aFloat! 

x2 

"return the number which is the current x location" 
rx2! 
x2: aFloat 

"save the number which is the current x location" 
x2 *- aFloat! 
xmax 

"return the number which is the maximum x location" 

txmax! 

xmax: aFloat 

"save the number which is the maximum x location" 

xmax «- aFloat! 
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xniin 

"return the number which is the minimum x location" 
txmin! 
xmin: aFloat 

"save the number which is the minimum x location" 
xmin «• aFloat! 

yi 

"return the number which is the previous y location" 
tyl! 
yl: aFloat 

"save the number which is the previous y location" 
yl «- aFloat! 

y2 

"return the number which is the current y location' 1 
ty2! 
yl: aFloat 

"save the number which is the currency location" 
y2 «■ aFloat! 
ymax 

"return the number which is the maximum y location" 

tyrnaxi 

ymax: aFloat 

"save the number which is the maximum y location" 

ymax «- aFloat! 

ymin 

"return the number which is the minimum y location" 

tymin! 

ymin: aFloat 

"save the number which is the minimum y location" 

ymin «- aFloat! 

zmax 

"return the number which is the maximum z location" 
tzmaxl 
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zmax: aFloat 

"save the number which is the maximum z location" 
zmax «- aFloat! 
zmin 

"return the number which is the minimum z location" 
tzmin! 
zmin: aFloat 

"save the number which is the minimum z location" 

zmin «• aFloat! ! 
IGerber methodsFon 'private*! 
execute 

"process the parsed Gerber command line" 



CM02' = self m) ifFalse: [ "not the end of file" 
CG54' a self g) ifTrue: [ "change the aperture" 
oldApt «- self d. 

apShape «- (apertures at: self d) shape, 
size *• (apertures at: self d) size,] 
ifFalse: [ "create the plot entity" 

('DOr = self d) ifTrue: [ "draw a trace" 

self makeTrace.]. 
CD02* = self d) ifTrue: [ " move to new location" 
self xl: self x2. 
self yl: selfyl]. 
CD03* = self d) ifTrue: [ "create a pad" 
self makePad.]. 

]. 

]! 

extents 

"return the absolute extents of the coordinates" 

I aPoint x y zl 

z «- zmax. 
x <-(xmax + 1.0). 
y *-(ymax + 1.0). 
taPoint «- x@y@z! 

makePad 

"create a pad and put it in the data base" 

I z aPoint aPad I 

aPoint «- x2@y2@0.0. 

aPad «- Pad new. 

aPad coordinates: aPoint, 

aPad padshape: self apShape. 

aPad padsize: self size. 
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self flashAdd: aPad. 
self xl:x2. 
self yl: y2I 

makeTrace 

"create a trace and put it in ihe data base" 

I z aPoint bPoint aTrace I 
bPoint«-x2@y2@0.0. 
aPoint <-xl@yl@0.0. 
aTrace «■ Trace new. 
aTrace beginPoint: aPoint. 
aTrace endPoint: bPoint. 
aTrace w: self size, 
self traceAdd: aTrace. 
selfxi:x2. 
self yl:y2! 

readFile 

"Read and process a Gerber file. 

This is the main method which kicks off the rest of the 

Gerber methods" 

I gerberFile commandLine aFileName aFile theFileNamel 

theFileName *■ FilllnTheBlank 

request:* What is the name of the aperture file?' 

initialAnswer: 'aperture.dict' . 
aFile «- FileStream oldFileNamed: theFileName. 

apertures *- (aFile fileln). "apertures becomes the dictionary containing all aperture information" 
aFile close. 

self xl: 0.0. "set all coordinates to zero" 

self yl: 0.0. 

self x2: 0.0. 

self y2: 0.0. 

self m: 'M00\ 

self xmin: 0.0. 

self xmax: 0.0. 

self ymin: 0.0. 

self ymax: 0.0. 

self zmin: 0.0. 

self zmax: 12.0. 

traces «- OrderedCollection new: 2. 
flashes «- OrderedCollection new:3. 

aFileName «- FilllnTheBlank 

request: 'What is the Gerber file name?* 

initialAnswer 'gecomp.gbr'. 
gerberFile *■ (FileStream oldFileNamed: aFileName) readonly. 
gerberFile reset. 
[gerberFile atEnd] whileFalse:[ 
commandLine «- String readFile: gerberFile. 
Transcript cr, show: commandLine; cn endEntry. 
self parse: commandLine. 
"self inspect." 
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self execute, 
"self inspect/ 
commandLine release.]. 
gerberFile close.! ! 
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Object subclass: # Aperture 

instance VariableNames: 'name shape size * 

class VariableNames: " 

poolDictionaries: " 

category: 'Graphics-Circuits'! 
Aperture comment: 

'This class is an aperture for a Gerber photoplotter'! 



! Aperture methodsFon 'comparing'! 

<= anAperture 

"compare the name of the apertures 
for an ASCII sort" 

t self name <= anAperture name! ! 

! Aperture methodsFon 'accessing'! 

name 

"return the string that is the Gerber D 
number for this aperture" 

t name! 

name: aString 

"save the string that is the Gerber D 
number for this aperture" 

name*- aString! 

shape 

"return the string that is the shape 
for this aperture" 

t shape! 

shape: aString 

"save the string that is the shape 
for this aperture" 

shape «■ aString! 

size 

"return the number that is the size 
for this aperture in inches" 

t size! 

size: aFloat 

"save the number that is the size 
for this aperture in inches" 

size <- aFloat! ! 
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Line3D subclass: #Trace 

instance VariableNames: 'w surface ' 

class VariabieNames: " 

pooIDictionaries: " 

category: ' Graphics-Circuits ' ! 
Trace comment: 
'This class is the structure for 
printed circuit traces. 
Traces are 3D lines with width. 
They are 3D because of the 
MCB application.'! 



!Trace methodsFor. 'accessing'! 
surface 

"Return the surface associated with the circuit trace" 

t surface! 

surface: aString 

"Set the surface associated with the circuit trace" 

t surface <-aString! 

w 

"Return the line width of the circuit trace" 
t w! 
w: value 

"Set the line width of the circuit trace" 

tw<- value! ! 

[Trace methodsFor. 'transformation'! 

convToView: a3DPoint 

"convert world coordinates to view coordinates based on world extents 
a3DPoint is the extents of the world view." 

I aTrace pointl point2 nTrace wPoint i 

wPoint- 1024.0@800.0@ 800.0. 

pointl «■ self beginPoint "begin translation from world to display coordinates" 
point2 *■ self endPoinL 

pointl «- (pointl /a3DPoint) * wPoint " constant is view maximums" 
point2 «- (point2 / a3DPoint) * wPoint. "assume z extent in view same as y" 
pointl y: ( 0.0 - (pointl y)). "invert y to move display origin to lower left" 
point2 y: ( 0.0 - (point2 y)). 
nTrace «■ Trace new. 
nTrace beginPoint: pointl. 
nTrace endPoint point2. 
"nTrace inspect" 
t nTrace! 
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mapToSphere: a3DPoint 

"map circuit to sphere centered in the model volume." 

I aTrace pointl point2 nTrace abcdtphxlylzlrl 

pointl «■ self beginPoint "begin translation from world to display coordinates" 

point2 *- self endPoint. 

a «- (Float pi). 

b«-(-1.0*(a/2.0)). 

c*-b. 

d*-(3.0*a)/4.0. 

a3DPoint z: 1.0. "we don't want to change z." 
r<-(a3DPoint x)/4.0. 

pointl *• (pointl / a3DPoint). " set coordinates to fractional part of volume." 

point2 (point2 / aSDPoint). 

xl «- pointl x. 

yl pointl y. 

zl «- pointl z, 

t«-(a*xl)+b. 

ph «- (c*yl)+<L 

xl «- (r*(t sin)*(ph sin))+((a3DPoint x)/2.0). 

yl «- (r*(ph cos))+((a3DPoint y)/2.0). 

zl «- (r*(tcos)*(ph sin))+((a3DPoint z)/2.0). 

zl<-zl + pointl z. 

pointl x: xl. 

pointl y: yl. 

pointl z: zl. 

xl «- point2 x. 

yl «- point2 y. 

zl «• point2 z. 

t«-(a*xl)+b. 

ph «• (c*yl)+<L 

xl «- (r*(t sin)*(ph sin))+ ((a3DPoint x)/2.0). 

yl «- (r*(ph cos))+((a3DPointy)/2.0). 

zl «■ (r*(t cos)*(ph sin))+((a3DPoint z)/2.0). 

zl «■ zl + point2 z. 

point2 x: xl. 

point2y: yl. 

point2 z: zl. 

pointl «■ pointl * a3DPoint. 
point2 «■ point2 * a3DPoint. 

pointl x: (pointl x +(a3DPoint x / 2.0)). "translate from origin back to center of volume" 
pointl y: (pointl y+(a3DPoint y / 2.0)). 

point2 x: (point2 x +(a3DPoint x / 2.0)). "translate from origin back to center of volume" 
point2 y: (point2 y+(a3DPoint y / 2.0)). 
self beginPoint pointl. 
self endPoint: point2. 
"self inspect"! 

rotate: aMatrix 

"rotate trace as defined by transformation matrix" 

I ml m2 1 

ml «- self beginPoint asMatrix. 
m2 «• self endPoint asMatrix. 
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ml «- ml * aMatrix. 

ml «• m2 * aMatrix. 

self beginPoint: ml asPoinGD. 

self endPoim: m2 asPoinGD! ! 
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DisplayObject subclass: #Pad 

instance VariableNames: 'padshape padsize coordinates surface * 

class VariableNames: " 

poolDiciionaries: " 

category: 'Graphics-Circuits'! 
Pad comment: 

'This class is the structure for 
constructing printed circuit pads. 
At this rime pads are restricted 
to the shapes round and square/! 



!Pad methodsFon 'accessing* I 

coordinates 

"Answer the center point on the receiver." 

tcoordinates! 

coordinates: a3DPoint 

"Answer the center point on the receiver." 

coordinates a3DPoint! 

padshape 

"return the string 'circle' or 'square'" 

tpadshape! 

padshape: aString 

"set to the string 'circle' or 'square'" 

padshape «- aString! 

padsize 

"return the size in inches, 
either a diameter or side." 

tpadsize! 

padsize: aNumber 

"set the size in inches, 
either a diameter or side." 

padsize*- aNumber! 

surface 

"Answer the surface associated with the receiver." 

t surface! 

surface: aString 

"Save the surface associated with the receiver." 

surface <- aString! ! 
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IPad methodsFon 'transformation'! 

convToView: a3DPoint 

"convert world coordinates to view coordinates based on world extents 
a3DPoint is the extents of the world view/* 

I aPad point 1 point2 nPad wPoint I 

wPoint- 1024.0@800.Q@ 800.0. 

pointl «- self coordinates, "begin translation from world to display coordinates" 
pointl «- (pointl / a3DPoint) * wPoint. " constant is view maximums" 

"assume z extent in view same as y" 
pointl y: ( 0.0 - (pointl y)). "invert y to move display origin to lower left" 
nPad <- Pad new. 
nPad coordinates: pointl. 
"nPad inspect." 
t nPad! 

mapToSphere: a3DPoint 

"convert world coordinates to view coordinates based on world extents 
a3DPoint is the extents of the world view." 

I aPad pointl point2 nPad a b c d t ph xl yl zl r I 

pointl *- self coordinates, "begin translation from world to display coordinates" 

a «■ (Roat pi). 

b<-(-1.0*(a/2.0)). 

c«~b. 

d<-(3.0*a)/4.0. 

a3DPoint z: 1.0. "we don't want to change z." 

pointl «- (pointl / a3DPoint). " constant is view maximums" 

"assume z extent in view same as y" 

r «- (a3DPoint x) /4.O. 
xl pointl x. 
yl «- pointl y. 
zl «■ pointl z. 
t«~(a*xl>fb. 
ph (c*yl)+± 

xl (r*(t sin)*(ph sin))+((a3DPoint x)/2.0). 

yl «- (r*(ph cos))+((a3DPoint y)/2.0). 

zl «- (r*(t cos)*(ph sin))+((a3DPoint z)/2.0). 

zl<-zl + pointl z. 

pointl x: xl. 

pointl y: yl. 

pointl z: zL 

self coordinates: pointl. 

"self inspect."! 

rotate: aMatrix 

"rotate pad as determined by transformation matrix" 

I ml I 

ml «- self coordinates asMatrix. 

ml «- ml * aMatrix. 

self coordinates: ml asPoint3D! ! 
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!Pad methodsFon 'initialize-release*! 

coordinates: a3DPoint padshape: aString padsize: aFloat 
"set the variables for a pad" 

self coordinates: a3DPoint* 
self padshape: aString. 
self padsize: aFloat! ! 
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Object subclass: #Iges 

instance VariableNames: 'diet param surfaces ' 

class VariableNames: " 

poolDiciionaries: " 

category: 'Graphics-Circuits'! 
Iges comment: 

'This class reads, writes, and processes IGES files/! 



Ilges methodsFon 'private*! 

readFile 

"Read and process an IGES file. 

This is the main method which kicks off the rest of the 

IGES methods" 

I igesFile commandLine commandLine2 aFileName aFile deNumber aString index I 

aFileName «• FiillnTheBlank 

request: 'What is the IGES file name?' 

initialAnswer: 'surfacciges'. 
igesFile *■ (FileStream oldFileNamed: aFileName) readonly. 
igesFile reset. 

[igesFile atEnd] whileFalse:[ 
commandLine «- String readRecord: igesFile. 

Transcript cr, show:'size ^print:(commandLine size);show: ' type chr ';print: (commandLine at: 
73);cr,endEntry. 

Transcript cr, show: commandLine; cr, endEntry. 
(commandLine at: 73) — SD 
ifTrue:[ 

commandLine2 *• String readRecord: igesFile. 
Transcript cr, show:'size ';print:(comrnandLine size);show: ' type chr ';print: (commandLine at: 
73);cr,endEntry, 

self getD: commandLine with: commandLine2] 
ifFalse: [(commandLine at: 73) = SP 
ifTrue:[ 

index <- 65. 

[(commandLine at: index) = S ] 
whileTrue: [index index + 1]. 

deNumber «- ((commandLine copyFrom: index to: 72) asNumber). 

self getP: commandLine with: deNumber]]. 
commandLine release.]. 
igesFile close.! ! 

Ilges methodsFon 'accessing'! 

diet 

"Return the dictionary which is the directory portion of the IGES file" 
t diet! 
dice aDictionary 

"Create the dictionary which is the directory portion of the IGES file" 
t diet «- aDictionary! 
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param 

"Return the collection which is the parameter geometry portion of the IGES file" 

t param! 

param: aCoilection 

"Create the collection which is the parameter portion of the IGES file" 

t param «- ^Collection! 

surfaces 

"Return the surfaces which are created from the IGES file" 

t surfaces! 

surfaces: aCoilection 

"Create the collection which will be the surfaces created from the IGES file" 

t surfaces «- aCoilection! ! 

!Iges methodsFon 'parsing*! 

getD: aString with: aString2 

"process the directory entry into dictionary keyed by directory number" 

I dir newString ! 
dir «• Directory new. 

dir eType: ((newString «- self getField: aString at: 1) asNumber). 

dir paramDataPntn ((newString «- self getField: aString at: 9) asNumber). 

dir refByOthen ((newString «• self getField: aString at: 65) asNumber). 

dir dirNumben ((newString self getField: aString at 74) asNumber), 

dir paramRecCnc ((newString «- self getField: aString2 at: 25) asNumber). 

dir eLabel: (newString «- self getField: aString2 at: 57). 

"dir inspect" 

self diet at: (dir dirNumber) put: dir! 

getField: aString at: aNumber 

"get the required field from the command line string" 

I stop newString index I 

index <- aNumber. 

(aString size - aNumber) < 7 

ifTrue: [stop <- (aString size - aNumber)] 

ifFalse: [stop «- 7]. 
[(aString ac index) = S ] 
whileTrue:[ index *- index + 1 ]. 
index > (aNumber + stop) 
ifTrue: [ newString *• '0'] 
ifFalse: [ 

newString «- aString copyFrom: index to: (aNumber + stop).]. 
"newString inspect." 
t newString! 

getLine: aString 
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"process a parameter command for a line" 

I xl x2 yl y2 zl il type stop index aLine newString I 

index «■ L 

stop «- aString indexOfSubCoIlection: Y startingAt: index. 

type «- (newString «• (aString copyFrom: index to: (stop-1)) asNumber). 

index «■ stop + 1, 

newString release. 

type= 110 

ifTrue: [ 

stop <- aString indexOfSubCoIlection: 7 startingAt: index. 

xl «■ (newString «- (aString copyFrom: index to: (stop-1)) asNumber). 

index «- stop + 1. 

newString release. 

stop «- aString indexOfSubCoIlection: Y startingAt: index. 

yl «• (newString «- (aString copyFrom: index to: (stop-1)) asNumber). 

index «• stop + 1. 

newString release. 

stop «- aString indexOfSubCoIlection: Y startingAt: index. 

zl «- (newString «■ (aString copyFrom: index to: (stop-1)) asNumber). 

index «- stop + 1. 

newString release. 

stop «- aString indexOfSubCoIlection: 7 startingAt: index. 

x2 «• (newString «- (aString copyFrom: index to: (stop-1)) asNumber). 

index «■ stop + 1. 

newString release. 

stop «- aString indexOfSubCoIlection: Y startingAt: index. 

y2 «- (newString «- (aString copyFrom: index to: (stop-1)) asNumber). 

index <- stop + 1. 

newString release. 

stop <- aString indexOfSubCoIlection: Y startingAt: index. 

22 «• (newString «- (aString copyFrom: index to: (stop-1)) asNumber). 

aLine +- Line3D new. 

aLine beginPoinc xl@yl@zl. 

aLine endPoint: x2@y2@z2. 

t aLine] 

ifFalse: [ Transcript cr t show: 'Not a Line'; cr, 

show: 'Entity type = *; print: type;cn endEntry.]! 

getP: aString with: dirNumber 

"process the parameter file for lines or surfaces" 

I index stop type newString aLine aSurface i 

index «- 1. 

stop «- aString indexOfSubCoIlection: Y startingAt index, 

type «• (newString «- (aString copyFrom: index to: (stop-1)) asNumber). 

index*- stop + 1. 

newString release. 

type= 110 

ifTrue: [ aLine «■ self getLine: aString. 

self param at: dirNumber put: aLine.] 
ifFalse: [aS urface «• self getSurface: aS tring. 

aSurface name: ((diet an dirNumber) eLabel). 

self surfaces add: aSurface.]! 
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getSurface: aString 

"process a parameter command for a surface" 

I el e2 stop index aSurface type newString anEdge anEdge2 I 
index *■ 1. 

anEdge <- Line3D new. 
anEdge2 «• Line3D new. 

stop aString indexOfSubCollection: 7 startingAt: index. 

type «■ (newString «• (aString copyFrom: index to: (stop-1)) asNumber). 

index <- stop + 1. 

newString release. 

type= 118 

ifTrue:[ 

stop <- aString indexOfSubCollection: 7 startingAt index. 

el «• (newString «- (aString copyFrom: index to: (stop-1)) asNumber). 

index* stop + 1. 

newString release. 

stop *• aString indexOfSubCollection: V startingAt: index. 
e2 «• (newString *• (aString copyFrom: index to: (stop-1)) asNumber). 
anEdge beginPoint: ((param at: el) beginPoint deepCopy). 
anEdge endPoint: ((param at: el) endPoint deepCopy). 
anEdge2 beginPoint: ((param at; e2) beginPoint deepCopy). 
anEdge2 endPoint: ((param at: e2) endPoint deepCopy). 
aSurface «■ Surface el: anEdge e2: an£dge2. 
t aSurface] 

i£False:[ Transcript cn show: 'Not a Surface'; en 

show: 'Entity type = '; print: type;cr; endEntry.]! ! 
„ ._ — ~ . _ „ "I 

Iges class 

instanceVariableNames: "! 



llges class methodsFon 'instance creation'! 
new 

"Set up IGES for holding graphic entities created as file is processed" 
! newSelf I 

newSelf*- super new. 

newSelf diet: Dictionary new. 

newSelf param: Dictionary new. 

newSelf surfaces: (SortedCollection sortBlock: 

[-.surface 1 :surface2l 

surface 1 name <= surface2 name]), 
t newSelf! ! 
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Object subclass: #Directory 

instanceVariableNames: 'eLabei eType paramDataPntr paramRecCnt dirNumber refByOther 1 

ciassVariableNames: " 

poolDictionaries: *' 

category: 'Graphics-Circuits'! 
Directory comment: 

'This class is used to hold the directory entry 
information from IGES files.' i 



[Directory methodsFon 'accessing'! 
dirNumber 

"return the First directory number for the graphic entity" 

t dirNumber! 

dirNumber aNumber 

"save the first directory number for the graphic entity" 

dirNumber «- aNumber! 

eLabei 

"return the name of the graphic entity" 

teLabel! 

eLabei: aString 

"save the name of the graphic entity" 

eLabei <-aString! 

eType 

"return the numeric type of the graphic entity" 

teType! 

eType: aNumber 

"save the numeric type of the graphic entity" 

eType «• aNumber! 

paramDataPntr 

"return the number of the first record of the graphic entity 
in the paramater section" 

TparamDataPntr! 

paramDataPntr aNumber 

"save the number of the first record of the graphic entity 
in the paramater section" 

paramDataPntr «• aNumber! 

paramRecCnt 
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"return the number of records for the graphic entity 
in the paramater section" 

tparamRecCnt! 

paramRecCnt:aNumber 

"save the number of records for the graphic entity 
in the paramater section" 

paramRecCnt *- aNumber! 

refByOther 

"return if the graphic entity referenced by other entities" 

trefByOther! 

refByOthenaNumber 

"save status if the graphic entity referenced by other entities" 

refByOther «- aNumber! ! 
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Object subclass: #Surface 

instance VariableNames: 'name el e2 e3 e4 matrix normal edgenorm ce sina cosa sinb cosb cirot r ri 

classVariableNames: " 
pooIDicuonaries: " 
category: 'Graphics-Circuits'! 
Surface comment: 

'Each surface is a coplanar ruled surface with straight lines 
as its edges. Each edge is a member of the Line3D class. Each 
end point on an edge is a vertex. Each vertex is a member of the 
PoinGD class/! 



!Surface methodsFon 'accessing'! 
ce 

"return common edge ce" 

t ce! 

ce: aString 

"save common edge ce" 

ce*- aString! 

cirot 

"return the circuit rotation matrix " 

t cirot! 

cirot: aMatrix 

"save the circuit rotation matrix" 

cirot <- aMatrix! 

cosa 

"return cos of angle" 

t cosa! 

cosa: aNumber 

"save cos of angle" 

cosa «• aNumber! 
cosb 

"return cos of angle" 

t cosb! 

cosb: aNumber 

"save cos of angle" 

cosb aNumber! 

el 
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"return edge 1" 

tel! 

el: aLine3D 

"store edgeT 

el «• aLine3D! 

e2 

"return edge 2" 

te2! 

e2: aLine3D 
"store edge 2" 

e2 «- aLine3D! 

e3 

"return edge 3" 

te3! 

e3: aLine3D 

"store edge 3" 

e3 *- aLine3D! 

e4 

"return edge 4" 

te4! 

e4: aLine3D 

"store edge 4" 

e4 *• aLine3D! 

edgenorm 

"return the array of inward edge normals" 

t edgenorm! 

edgenorm: an Array 

"save the array of inward edge normals" 

edgenorm <- an Array! 

matrix 

"return transformation matrix" 

t matrix! 

matrix: aMatrix 

"store transformation matrix" 
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matrix <- aMatrix! 

name 

"return name" 

t name! 

name: aString 

"save name of surface" 

name<- aString! 

normal 

"return normal vector of surface" 

t normal! 

normal: aVector 

"save surface normal vector" 

normal «- aVector! 

r 

"return the surface rotation matrix" 

tr! 

r: aMatrix 

"save the surface rotation matrix" 

r «- aMatrix! 

ri 

"return the surface inverse rotation matrix" 

tri! 

ri: aMatrix 

"save the surface inverse rotation matrix" 

ri<- aMatrix! 

sina 

"return sin of angle" 

t sina! 

sina: aNumber 

"save sin of angle" 

sina «- aNumber! 

sinb 

"return sin of angle" 
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t sinb! 

sinb: aNumber 

"save sin of angle" 

sinb «- aNumber! ! 

ISurface methodsFon 'testing'! 

commonEdge: aSurface 

"fine the common edge with another surface" 

(self el sameAs: aSurface el) ifTrue:[ self ce: *el\]. 
(self el sameAs: aSurface e2) ifTnie:[ self ce: 'el'.], 
(self el sameAs: aSurface e3) ifTrue:[ self ce: 'el'.], 
(self el sameAs: aSurface e4) ifTrue:[ self ce: 'el'.], 
(self e2 sameAs: aSurface el) ifTrue:[ self ce: 'e2\]. 
(self e2 sameAs: aSurface e2) ifTrue:[ self ce: *e2\]. 
(self e2 sameAs: aSurface e3) ifTrue:[ self ce: *e2\], 
(self e2 sameAs: aSurface e4) ifTrue:( self ce: *e2\]. 
(self e3 sameAs: aSurface el) ifTrue:[ self ce: 'e3'J. 
(self e3 sameAs: aSurface e2) ifTrue:[ self ce: 'e3'.]. 
(self e3 sameAs: aSurface e3) ifTrue:[ self ce: , e3 , J. 
(self e3 sameAs: aSurface e4) ifTrue:[ self ce: *e3'.]. 
(self e4 sameAs: aSurface el) ifTrue:[ self ce: 'e4\]. 
(self e4 sameAs: aSurface e2) ifTrue:[ self ce: 'e4\]. 
(self e4 sameAs: aSurface e3) ifTrue:[ self ce: f e4' J. 
(self e4 sameAs: aSurface e4) ifTrue:[ self ce: f e4'.]! ! 

[Surface methodsFon 'private*! 

clipPad: aCollectionOfPads 

"Clip all pads into two lists. Those that fit on the surface 

and those parts which don't. Cyrus-Beck 2D clipping algorithm." 

I internalCollection externalCollection cnt w wnp col aPad I 
internalColiection «- OrderedCollection new: L 
externalCollection *■ OrderedCollection new: 1. 
w Array new: 4. 
col ** Array new: 2. 

1 to: (aCollecuonOfPads size) do:[ :index i 
aPad*- aCollecuonOfPads at: index, 
cnt «■ 0. 

w at: 1 put: (aPad coordinates - self el beginPoint) as Vector, 
w at: 2 put: (aPad coordinates - self e2 endPoint) as Vector, 
w au 3 put: (aPad coordinates - self el beginPoint) as Vector, 
w at: 4 put: (aPad coordinates - self e2 endPoint) as Vector. 
1 to: 4 do: [:edg I 

wnp «- (w at: edg) dotProduct: (self edgenorm at: edg). 
"Transcript cr,show:Tad \print: index;show:' w ^printtfw at: edg);show:' wnp *;prim: wnp;cnendEntry. M 
(wnp < 0) 

ifTrue:[cnt«-cnt + I.]. 
]- 

(cnt = 0) ifTrue:[ aPad surface: self name. 

internalCollecdon add: aPai] 
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ifFalse: [externalCollection add: aPad.]. 

]. 

col at: 1 put: internalCollection. 
col at: 2 put: externalCollection. 
t col! 

clipTrace: aCollectionOfTraccs 

"Clip all traces into two lists. Those that fit on the surface 

and those parts which don't. Cyrus-Beck 2D clipping algorithm." 

I internalCollection externalCollection tl tu i k d w 
aTrace aTracel aTrace2 aTrace3 dnp wnp tpl p2 col i 
internalCollection «• OrderedCoIlecdon new: 1. 
externalCollection *• OrderedCoilection new: L 
w «- Array new: 4. 
col «- Array new: 2. 
aTracel <- Trace new. 
aTrace2 *■ Trace new. 
aTrace3 «- Trace new. 

1 to: (aCollectionOfTraces size) do:[ :index I 
aTrace aCollectionOflraces at: index. 
tl<-0. 
tu<- 1. 
k<-0. 

d <- aTrace asVector. 

w at: 1 put: (aTrace beginPoint - self el beginPoint) asVector. 
w at: 2 put: (aTrace beginPoint - self e2 beginPoint) asVector. 
w at: 3 put: (aTrace beginPoint - self e3 beginPoint) asVector. 
w at: 4 put: (aTrace beginPoint - self e4 beginPoint) asVector. 
"self halt." 

1 to: 4 do: [:edg I 

dnp <- d dotProduct: (self edgenorm at: edg). 

wnp «- (w at: edg) dotProduct (self edgenorm at: edg). 

(dnp = 0.0) ifFalse:[ 

t«- -1.0*(wnp/dnp). 

"Transcript cnshow:* dnp Sprint dnp;show:' wnp ';print: wnp; show:' t *;print: t;cr t endEntry.' 
"Transcript cr;show:'d ^printd^how:' w *;print:(w an edg);show:' enorm ';print: (self 
edgenorm at: edg);cr;endEntry." 
(dnp > 0) 
ifTrue:[(t> 1) 

ifFalse:[ tl*- d max: t] 
ifTrue:[k*-l.].] 
ifFalse:[(t<0) 

i£False:[ tu «- tu min: t] 
ifTrue:[k«-L].]. 

] 

i£True:[ (wnp < 0) ifTrue:[k «• 1.].].]. 
"Transcript cr t show: 'index ';print: index ; cr t 

show:'tl sprint: tl; cr, 

show: 'tu \print tu;cnendEntry." 
(k=l)ifFaise:[ 
(tl>tu) 

ifFalse:[((d=0)&(tu^l)) "if true whole trace is on surface" 
ifTrue:[ aTrace surface: self name. 

" set trace as surface dependent?" 
internalCollection add* aTrace deepCopy.] 
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ifFaise:[ 

pi «- ((aTrace beginPoint)+(((aTrace endPoint)-(aTrace beginPoint))*tl)). 
p2 «- ((aTrace beginPoint)+(((aTrace endPoint)-(aTrace beginPoint))* tu)). 
((tl > 0) & (tu < 1)) rt if true split line into 3 segments" 
ifTrue: [aTracel beginPoint: aTrace beginPoint. 

aTracel endPoint: pi. 

aTracel w: aTrace w, 

aTrace2 beginPoint: pi. 

aTrace2 endPoint: p2. 

aTrace2 w: aTrace w. - 

aTrace2 surface: self name. 

aTrace3 beginPoint: p2. 

aTrace3 endPoint: aTrace endPoint. 

aTrace3 w: aTrace w. 

internalCollection add: aTrace2 deepCopy. 
externalCollection add: aTracel deepCopy. 
externalCollection ad± aTrace3 deepCopy.] 
ifFalse:[((U=0)&(tu<l)) 

ifTrue:[aTrace2 beginPoint: aTrace beginPoint. 

aTrace2 endPoint: p2. 

aTrace2 w: aTrace w. 

aTrace2 surface: self name. 

aTrace3 beginPoint: p2. 

aTrace3 endPoint: aTrace endPoint. 

aTrace3 w: aTrace w. 

internalCollection add: aTrace2 deepCopy. 

externalCollection ad± aTrace3 deepCopy.] 
ifFalse:[aTrace2 beginPoint: pL 
aTrace2 endPoint: aTrace endPoint. 
aTrace2 w: aTrace w. 
aTrace2 surface: self name. 
aTrace3 beginPoint: aTrace beginPoint 
aTrace3 endPoint: pi. 
aTrace3 w: aTrace w. 
internalCollection add: aTrace2 deepCopy. 
externalCollection ad± aTrace3 deepCopy.],]. 
].]. 

] 

ifTme: [externalCollection add: aTrace deepCopy.]. 

]. 

col an 1 put: internalCollection. 
col at: 2 put: externalCollection. 
t col! 



flndAngle: aSurface 

"find angles for rotation between surfaces" 



I vl v2 v3 v4 x yi 

x <- Vector x: L0 y: 0.0 z: 0.0 w: 0.0. 
y «- Vector x: 0.0 y: 1.0 z: 0.0 w: 0.0. 
vl self normal cos: aSurface normal. 
v2 «- self normal sin: aSurface normal. 
v3 «■ x cos: aSurface normal. 
v4 <- x sin: aSurface normal. 

"Transcript cr;show:'sina Sprint: v2;show:' cosa ^print:vl;cr,endEntry. 
Transcript cr,show:'sinb ';print: v4;show:' cosb T ;print:v3;cr;endEntry. M 
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self cosa: vl. 
self sina: v2. 
self cosb: v3. 
self sinb: v4! 

findNormal 

"calculate the inward normal vector of an edge" 

I vl v2 v3 v4 pi p2 p3 p4 ql q2 q3 q4 anArrayl 

anArray *- Array new: 4. 
vl ♦-self el as Vector. 
v2 «• self e2 asVector. 
v3 «- self e3 asVector. 
v4 *- self e4 asVector. 

pi *- (PoinGD x:(-1.0*vl y) y: (vl x) z:(vl z)) asVector. 
p2 «- (PoinGD x:(-1.0*v2 y) y: (v2 x) z:(v2 z)) asVector. 
p3 «- (PoinGD x:(-1.0*v3 y) y: (v3 x) z:(v3 z)) asVector. 
p4 #- (PoinGD x:(-1.0*v4 y) y: (v4 x) z:(v4 z)) asVector. 
ql «- ((self e2 endPoint) - (self el beginPoint)) asVector. 
q2 *• ((self el endPoint) - (self e2 beginPoint)) asVector. 
q3 *■ ((self e2 endPoint) - (self e3 beginPoint)) asVector. 
q4 ((self e2 beginPoint) - '(self e4 beginPoint)) asVector. 
((pi dotProduct: ql) < 0) 
ifTrue:[pl x: (-1.0 *(pl x». 

pi y: (-1.0*(pl y)).]. 
((p2 dotProduct: q2) < 0) 
ifTrue:[p2 x: (-L0 *(p2 x)). 

p2y: (-L0*(p2 y)).]. 
((p3 dotProduct: q3) < 0) 
ifTrue:[p3 x: (-1.0 *(p3 x)). 

p3 y: (-L0*(p3 y))J. 
((p4 dotProduct: q4) < 0) 
ifTrue:[p4 x: (-1.0 *(p4 x)). 

p4y: (-L0*Gp4y))J. 
anArray at: 1 put (pi unit). 
anArray at: 2 put: (p2 unit). 
anArray at: 3 put: (p3 unit). 
anArray at: 4 put: (p4 unit), 
r anArray! 

form: aForm 

"set the form for line width for the 2D edge lines" 

self el:(self el form: aForm). 
self e2:(self e2 form: aForm). 
self e3:(self e3 form: aForm). 
self e4:(self e4 form: aForm). 
tselfl! 

[Surface methodsFor: 'transforming'! 

buildRotX: costh with: sinth 

"build a rotation matrix for the X axis." 

Iranal 
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a «• Array new: 2. 

m «- Matrix new:4 by: 4. 

m sctToZero. 

m atPoint: 1@1 put: 1.0. 

m atPoint: 2@2 put: costh. 

m atPoint: 3@2 put: sinth. 

m atPoint: 2@3 put:(-1.0 * sinth). 

rn atPoint: 3@3 put: costh. 

m atPoint: 4@4 put: 1.0. 

n «- m fromPoint: 1@1 toPoinc 4@4. 

n atPoint: 3@2 put:(-1.0 * sinth). 

n atPoint: 2@3 put: sinth. 

a at: 1 put: m. 

a at: 2 put: n. 

ta! 

buildRotY: costh with: sinth 

"build a rotation matrix for the Y axis." 

Imnal 

a «- Array new: 2. 

m «- Matrix new:4 by: 4. 

m setToZero. 

m atPoint: 1@1 put: costh. 

m atPoint: 3@1 put: (-1.0 * sinth). 

m atPoint: 2@2 put: 1.0. 

m atPoint: 1@3 put:sinth. 

m atPoint: 3@3 put: costh. 

m atPoint: 4@4 put: 1.0. 

n «- m fromPoint: 1@ 1 toPoint: 4@4. 

n atPoint: 3@1 put: sinth. 

n atPoint: 1@3 put:(-1.0 * sinth). 

a at: 1 put m. 

a at: 2 put: n. 

ta! 

buildRotZ: costh with: sinth 

"build a rotation matrix for the Z axis." 

Imnal 

a «- Array new: 2. 

m «- Matrix new:4 by: 4. 

m setToZero. 

m atPoint: 1@ 1 put: costh. 

m atPoint: 2@1 put: sinth. 

m atPoint: 1@2 put: (-1.0 * sinth). 

m atPoint: 2@2 putxosth. 

tn atPoint: 3@3 put: 1.0. 

m atPoint: 4@4 put: 1.0. 

n «- m fromPoint: 1@ 1 toPoinn 4@4. 

n atPoint: 2@ 1 put: (-1.0 * sinth). 

n atPoinc 1@2 put: sinth. 

a at: 1 put: m. 
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a at: 2 put: n. 
t a! 

buiidRotZc: aMatrix 

"build a special circuit rotation matrix for the Z axis." 

I m sen a i 

a «■ Array new: 2. 

m * aMatrix fromPoint 1@ 1 toPoint 4@4. 

s «- m atPoint 2@ 1 deepCopy. 

c «- m atPoint 1@1 deepCopy. 

m atPoint: 1@1 put: s. 

m atPoint: 2@2 put: s. 

m atPoint: 2@ 1 put: c. 

m atPoint: 1@2 put: (-1.0 * c). 

n «- m fromPoint 1@ 1 toPoint 4<2>4. 

n atPoint: 2@ 1 put: ((n atPoinu 2@1)*-1.0). 

n atPoint: l@2put ((n atPoint: 1@2)*-1.0). 

a at: 1 put m. 

a at: 2 put: n. 

ta! 

buildTrans: aPoint3D 

"build a translation matrix from a 3D point" 

ItmaJ 

a «■ Array new: 2. 

t «- Matrix new:4 by: 4. 

t setToZero. 

t atPoint 1@4 put: ((aPoinGD x) * -1.0). 

t atPoint 2@4 put ((aPoint3D y) * -L0). 

t atPoint: 3@4 put (( aPoint3D z) * -1.0). 

t atPoint 4@4 put: 1.0. 

t atPoint: 1@1 put 1.0. 

t atPoint 2@2 put: 1.0. 

tatPoint 3@3 put 1.0. 

m «- 1 fromPoint 1@1 toPoint: 4@4. 

m atPoint 1@4 put(aPoint3D x). 

m atPoint 2@4 put:(aPoint3D y). 

m atPoint 3@4 put:(aPoint3D z). 

a at 1 put t 

a at: 2 put m. 

tal 

calcRotation: aSurface 

"calculate circuit and surface rotations for mapping to aSurface" 

I I ti tc tci rx rxi ry ryi rz rzi rzc rzci rztmp tmp I 
tmp *■ self buildTrans: (aSurface el beginPoint). 
t «• tmp at 1. 

ti «- tmp at 2. 

((self ce- •elO I (self ce='e2')) 
ifTrue:[ 

(selfce^'elOifTruezf 
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tmp self buUdTrans:(self el beginPoint).] 
i£False:[ tmp «- self buildTrans:(self e2 beginPoint).]. 
tc*- tmp at: L 
tci «- tmp at: 2. 

tmp <- self buildRotX: self cosa with: self sina. 
rx *■ tmp at: 1. 
rxi «- tmp at: 2. 
"self halt." 

(self cosb<0)ifTrue:[ 

self sinb: self sinb * -1.0. 

self cosb: self cosb * -1.0.]. 
((self sinb « 1) I (self sinb - -1)) ifTrue:[ self cosb: 1.0. 

self sinb: 0.0.]. 
tmp «■ self buildRotZ: self cosb with: self sinb. 
rz «- tmp at: 1. 
rzi «- tmp at: 2. 

( (aSurface normal y) >= 0) iflrue:[ 

self cirot: tc * rz * rxi * tci. " the order is crucial !!!!!!" 

self n t * rz * rx. 

self ri: rxi * rzi * ti.] 
ifFalse:[ 

self cirot: tc * rzi * rx * tci. 

self n t * rzi * rxi. 

self ri: rx * rz * ti.]. 
"self halt"] 
i£Faise:[ 

(self ce « , e3 f ) ifTrue:[ 

tmp «- self buildTrans:(self e3 beginPoint).] 
ifFalse:[ tmp «- self buildTrans: (self e4 beginPoint).]. 
tc«- tmp at: 1. 
tci <- tmp at: 2. 

tmp «- self buildRotY: self cosa with: self sina. 

ry ♦* tmp at: L 

ryi *• tmp at: 2. 

(self sinb<0) ifTrue:[ 

self sinb: self sinb * -1.0. 

self cosb: self cosb * -1.0.]. 
(self cosb = -1) ifTrue:[ self cosb: self cosb * -1.0.]. 
tmp «- self buildRotZ: self cosb with: self sinb. 
rz <- tmp at: 1. 
rzi «- tmp at: 2. 

( (aSurface normal x) >= 0) ifTnie:[ 

self cirot: tc * rzi * ry * tci. "the order is crucial !!!!!!" 

selfr t* rzi* ryi. 

self ri: ry * rz * ti.] 
ifFalse:[ 

self cirot: tc * rzi * ryi * tci. 

self n t * rz * ry. 

self ri: ryi* rzi* ti.]. 
"self halt."].! 

convToView: a3DPoint 

"convert world coordinates to view coordinates based on world extents 
a3DPoint is the extents of the world view." 

I aSurface pointl point2 nTrace wPoint aLine aLine2l 
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aLine «- Line3D new. 

aLine2 «• Line3D new. 

wPoint <- 1024.0@800.0@800.0. 

pointl *- self el beginPoint. "begin translation from world to display coordinates" 
point2 «- self el endPoint. 

pointl «- (pointl / a3DPoint) * wPoint. " constant is view maximums" 
point2 *• (point2 / a3DPoint) * wPoint. "assume z extent in view same as y" 
pointl y; ( 0.0 - (pointl y)). "invert y to move display origin to lower left" 
point2 y: (0.0 - (point2 y)). 
aLine beginPoint: pointl. 
aLine endPoint point2. 

pointl «- self e2 beginPoint. "begin translation from world to display coordinates" 
poim2 «• self e2 endPoint. 

pointl <- (pointl / a3DPoint) * wPoint. " constant is view maximums" 

point2 «- (point2 / a3DPoint) * wPoint. "assume z extent in view same as y" 

pointl y: ( 0.0 - (pointl y)). "invert y to move display origin to lower left" 

point2 y: (0.0 - (point2 y)). 

aLine2 beginPoint: pointl. 

aLine2 endPoint point2. 

aSurface «■ Surface el: aLine e2: aLine2. 

t aSurface! 

mapToNextSurface: traceList with: padList and: surfaces and: surflist 
"This is the heart of mapping program. Clip all circuitry 
for current circuitry. The according to surflist, map remaining circuitry 
to attached surfaces. Include their mapped circuitry in current list. 
Return mapped circuitry to calling object" 

I mappedCircuits internalTraces extemalTraces internalPads 
extemalPads tempCircuit tempTrace tempPad linkList nextSurface sname I 
mappedCircuits «- Array new:2. 

Transcript en show: 'Mapping surface: Sprint: self name;cr;endEntry. 

self edgenorm: ( self findNormal). 

tempCircuit «- self clipTrace: traceList 

internalTraces «• (tempCircuit at l)deepCopy. 

extemalTraces <- (tempCircuit at 2)deepCopy. 

tempCircuit «- self clipPad: padList 

internalPads «- (tempCircuit at: l)deepCopy. 

extemalPads «- (tempCircuit at: 2)deepCopy. 

"self halt" 

1 to:(surflist size) do:[:index 1 
((surflist at index) = nil) ifFaise:[ 
(self name = ((surflist at index) first name)) i£True:[ 
linkList «■ surflist at: index. 
2 to: linkList size do:[:cnt 1 

sname «- (linkList at cnt) name. 
1 to: (surfaces size) do:[:scnt i 
((surfaces at sent) name = sname) ifTnic:[ 
nextSurface «- surfaces at sent] J. 
self findAngle: nextSurface. 
self commonEdge: nextSurface. 
self calcRotation: nextSurface. 
tempTrace «- extemalTraces deepCopy. 
1 to: (tempTrace size) do:[:tl 

(tempTrace at t) rotate: self cirot]. 
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tempPad «- extemaiPads deepCopy. 
1 to: (tempPad size) do:(:ti 

(tempPad at: t) rotate: self cirot]. 
1 to: (surfaces size) do:[:ti 
(surfaces at: t) rotate: self rj. 
1 to: (tempTrace size) do:[:tl 

(tempTrace at: t) rotate: self rj. 
1 to: (tempPad size) do:[:tl 

(tempPad at: t) rotate: self r.]. 
tempCircuit «• nextSurface mapToNextSurface: tempTrace 
with: tempPad and: surfaces and: surflisL 
1 to: (surfaces size) do:[:tl 
(surfaces an t) rotate: self ri.]. 
tempTrace «- tempCircuit at: 1. 
tempPad <- tempCircuit at: 2. 
1 to: (tempTrace size) do:[:ti 

(tempTrace ac t) rotate: self ri.]. 
1 to: (tempPad size) do:[:tl 

(tempPad at: t) rotate: self ri.]. 
internalTraces addAllLast: tempTrace deepCopy. 
internalPads addAllLast: tempPad deepCopy. 

]- 

]JJ. 

mappedCircuits at: 1 put: internalTraces deepCopy. 
mappedCircuits at: 2 put: internalPads deepCopy. 
t mappedCircuits! 

rotate: aMatrix 

"rotate a surface as defined by a transformation matrix" 

I ml m2 m3 m4 m5 I 
ml <- self el beginPoint asMatrix. 
m2 *- self el endPoint asMatrix. 
m3 «- self e2 beginPoint asMatrix. 
m4 «- self e2 endPoint asMatrix. 
m5 «- self normal asMatrix. 
ml <-ml*aMatrix. 
m2 *~ m2*aMatrix. 
m3«- m3*aMatrix. 
m4«- m4* aMatrix. 
m5 «- m5 * aMatrix. 
self el beginPoint: ml asPoint3D. 
self e3 beginPoint: ml asPoint3D. 
self el endPoint: m2 asPoint3D. 
self e4 beginPoint: m2 asPoint3D. 
self e2 beginPoint: m3 asPoint3D. 
self e3 endPoint: m3 asPoint3D. 
self e2 endPoint: m4 asPoint3D. 
self e4 endPoint: m4 asPoint3D. 
self normal: m5 as Vector! 

scaleBy: aPoint 

"scale 2D edges of a surface for viewing" 

self el: (self el scaleBy: aPoint). 
self e2: (self e2 scaleBy: aPoint). 
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self e3: (self e3 scaleBy: aPoint). 
self e4: (self e4 scaleBy: aPoint). 
t self! ! 

[Surface methodsFon 'displaying'] 

displayOn: aDisplayMedium at: aPoint clippingBox: clipRect 

"The form associated with this Surface will be displayed, according 
to one of the sixteen functions of two logical variables (rule), at 
each point on the Line. Also the source form will be first ANDed 
with aForm as a mask. Does not effect the state of the Path." 



self el displayOn: aDisplayMedium at: aPoint clippingBox: clipRect 
self e2 displayOn: aDisplayMedium at: aPoint clippingBox: clipRect 
self e3 displayOn: aDisplayMedium at: aPoint clippingBox: clipRect 
self e4 displayOn: aDisplayMedium at: aPoint clippingBox: clipRect! 

displayOn: aDisplayMedium at: aPoint clippingBox: clipRect rule: anlnteger mask: aForm 
"The form associated with this Surface will be displayed, according 
to one of the sixteen functions of two logical variables (rule), at 
each point on the Line. Also the source form will be first ANDed 
with aForm as a mask. Does not effect the state of the Path." 



self el displayOn: aDisplayMedium at: aPoint clippingBox: clipRect rule: anlnteger mask: aForm. 
self e2 displayOn: aDisplayMedium at: aPoint clippingBox: clipRect rule: anlnteger mask: aForm. 
self e3 displayOn: aDisplayMedium at: aPoint clippingBox: clipRect rule: anlnteger mask: aForm. 
self e4 displayOn: aDisplayMedium at: aPoint clippingBox: clipRect rule: anlnteger mask: aForm 

aDisplayMedium 

drawLine: el form 

from: el beginPoint + aPoint 

to: el endPoint + aPoint 

clippingBox: clipRect 

rule: anlnteger 

mask: aForm. 
aDisplayMedium 

drawLine: e2 form 

from: e2 beginPoint + aPoint 

to: e2 endPoint + aPoint 

clippingBox: clipRect 

rule: anlnteger 

mask: aForm. 
aDisplayMedium 

drawLine: e3 form 

from: e3 beginPoint + aPoint 

to: e3 endPoint + aPoint 

clippingBox: clipRect 

rule: anlnteger 

mask: aForm. 
aDisplayMedium 

drawLine: e4 form 

from: ©4 beginPoint + aPoint 

to: e4 endPoint + aPoint 

clippingBox: clipRect 
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rule: anlnteger 
mask: aForm"! ! 

Surface class 

instance VariableNames: "! 



! Surface class methodsFon 'instance creation'! 

el: aLine3D e2: aLine3D2 

"this takes the 2 edges of a ruled surface and constructs the rest of 
the information needed for a surface" 

I aSurface newEdge newEdge2 vl v2 mat I 

aSurface «- Surface new. 

aSurface el: aLine3D. 

aSurface e2: aLine3D2. 

newEdge «- Line3D new. 

newEdge2 «- Line3D new. 

newEdge beginPoint: aLine3D beginPoint. 

newEdge2 beginPoint: aLine3D endPoint. 

newEdge endPoint: aLine3D2 beginPoint. 

newEdge2 endPoint: aLine3D2 endPoint. 

aSurface e3: newEdge. 

aSurface e4: newEdge2. 

vl «• ((newEdge endPoint)-(newEdge beginPoint)) as Vector. 
v2 *■ ((aLine3D endPoint)-(aLine3D beginPoint)) as Vector. 

aSurface normal: ((v2 xProduct: vl) unit). tt right hand rule - positive z towards screen" 
mat *- Matrix new: 4 by: 4. 
mat setToZero. 

mat atPoint: 1 @ 1 put: (aSurface normal x), 
mat atPoint: 2 @ 2 put: (aSurface normal y). 
mat atPoint: 3 @ 3 put: (aSurface normal z). 
mat atPoint: 4 @ 4 put: (aSurface normal w). 
aSurface matrix: mat. 
t aSurface! ! 

! Surface class methodsFor: *new T ! 

new 

"create an instance of a surface' 1 
I aSurface I 

aSurface «- super new. 
t aSurface! ! 
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Link subclass: #SurfList 

instanceVariableNames: 'name ' 

classVariableNames: " 

poolDictionaries: " 

category: 'Graphics-Circuits'! 
SurfList comment: 

'This class is the elements of the surface list 
which orders the mapping process.*! 



ISurfList methodsFor 'printing'! 

printOn: aStream 

"print surface name" 

aStream nextPutAll: 'Surface name: '.name priniString! ! 
ISurfListmethodsFor: 'accessing'! 
name 

"retrieve the name of a surface" 

t name! 

name: aString 

"save the name of a surface" 

t name*- aString! ! 
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'From Smalltalk-80, Version 23 of 13 June 1988 on 21 August 1989 at 9:03:04 pm 1 ! 

Point subclass: #Point3D 

instance VariableNames: 'z 7 

class VariabieNames: 'RandomSource * 

pooiDictionaries: " 

category: 'Graphics-Primitives'! 



!Point3D methodsFon 'accessing'! 



tz! 

z: value 

t z *- value! ! 

IPoinGD methodsFon 'arithmetic'! 

* scale 

"Answer a new Point that is the product of the receiver and scale (which is a Point 
or Number).' 1 

I scalePoint I 

(scale isMemberOf: Point3D) 

i£False:[ scalePoint «■ scale@scale@ scale.] 

ifTrue:[ scalePoint*- scale as3DPoint]. 

tx * scalePoint x @ (y * scalePoint y) @ ( z * scalePoint z)! 

+ delta 

"Answer a new Point that is the sum of the receiver and delta (which is a Point 
or Number). " 

I deltaPoint I 

deltaPoint delta as3DPoint. 

rx + deltaPoint x @ (y + deltaPoint y) @ ( z + deltaPoint z)! 
- delta 

"Answer a new Point that is the difference of the receiver and delta (which is a Point 
or Number)." 

I deltaPoint I 

deltaPoint «- delta as3DPoint. 

tx - deltaPoint x @ (y - deltaPoint y) @ ( z - deltaPoint z)! 
/scale 

"Answer a new Point that is the quotient of the receiver and scale (which is a Point 
or Number)." 

I scalePoint I 

scalePoint «- scale as3DPoinL 

tx / scalePoint x @ (y I scalePoint y) @ ( z / scalePoint z)! 
//scale 

"Answer a new Point that is the quotient of the receiver and scale (which is a Point 
or Number)." 



56 



EP 0 416 568 A2 



I scaiePoint I 

scaiePoint «- scale as3DPoint 

tx // scaiePoint x @ (y II scaiePoint y) @ ( z // scaiePoint z)! 

abs 

"Answer a new Point whose x and y are the absolute values of the receiver's 
x and v." 

tPoinGD x: x abs y: y abs z: z abs! ! 
!Point3D methodsFor: 'converting'! 
as2DPoint 

"Answer with 2D version of self' 

I newPoint nx nyl 
nx*-selfx. 
ny «- self y. 
newPoint «- x@y. 
t newPoint! 

as3DPoint 

"Answer the receiver itself." 

tseif! 
asMatrix 

"Answer with a matrix version of selF 
I mat I 

mat <- Matrix new: 1 by:4. 
matatPoint: 1(2)1 put: self x. 
mat atPoint: 1@2 put: self y. 
mat atPoint: 1(2)3 put: self z. 
matatPoint: l@4put: 1.0. 
t mat! 

asVector 

"Answer with a vector version of selT 

I new Vector nx ny nzl 
nx«-self x. 
ny <-self y. 
nz self z. 

new Vector <- Vector x: nx y: ny z: nz w: 1.0. 
t new Vector! 

middle: otherPoint displaced: range 

"Answer a point in the middle, randomly displaced. 
Keep coordinates integral for speed only." 



57 



EP 0 416 568 A2 



I newX newY newZ I 
newX «- (x + oiherPoint x) // 2. 
new Y «- (y + otherPoint y) // 2. 
newZ «- (z + otherPoint z) // 2. 

newZ «- newZ + (RandomSource next - 0.5 * range) truncated, 
r Point3D x: newX y: newY z: newZ! ! 

!Point3D methodsFon 'printing'! 

priniOn: aStream 

"The receiver prints on aStream in terms of infix notation," 

xprintOn: aStream. 
aStream nextPut S@. 
y printOn: aStream. 
aStream nextPut S@. 
z printOn: aStream! ! 

!Point3D methodsFor: 'private'! 

setX: xValue setY: y Value setZ: zValue 
x «- x Value, 
y «- y Value, 
z *• zValuei ! 

!Point3D methodsFor 'comparing'! 

< aPoint 

"Answer whether the receiver is 'above and to the left and behind' of the 
argument, aPoint." 

t(x < aPoint x and: [y < aPoint y]) and: [z < aPoint z]! 
<= aPoint 

"Answer whether the receiver is 'neither below and to the right and in front' of the 
argument, aPoint." 

t(x <= aPoint x and: [y <= aPoint y]) and: [z <= aPoint z]! 
= aPoint 

"Answer whether the receiver is 'equal' to the 
argument, aPoint." 

t(x = aPoint x and: [y = aPoint y]) and: [z = aPoint z]! 
> aPoint 

"Answer whether the receiver is 'below and to the right and in front' of the 
argument, aPoinL" 

t(x > aPoint x and: [y > aPoint y]) and: [z > aPoint z]I 
>= aPoint 

"Answer whether the receiver is 'neither above and to the left and behind* of the 
argument, aPoint." 

t(x >= aPoint x and: [y >= aPoint y]) and: [z >= aPoint z]! ! 
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!Point3D methodsFon 'vector'! 
xpr± aPoini3D 

"Treat two 3D points as vectors and return their cross product 

as a 3D point. This will be the normal vector for the first two vectors." 

lijkxlylzll! ! 

!Point3D methodsFon 'transformation'! 

convToView: a3DPoint 

"convert world coordinates to view coordinates based on world extents 
a3DPoint is the extents of the world view." 

} pointl I 

pointl *■ (self / a3DPoint) * 1000@ 800(2)800. " constant is view maximums" 
pointl y: (800 - (pointl y)). "invert y to move display origin to lower left" 
t pointl I ! 

!Point3D methodsFor: 'copying'! 

deepCopy 

"Answer a copy of the receiver with its own copy of each instance variable." 

"Implemented here for better performance." 

tx deepCopy @ y deepCopy <§> z deepCopy! 

shallowCopy 

"Implemented here for better performance." 

tx@ y @ z! ! 

^ ^ _^ W| 



Point3D class 

instance VariableNames: "! 



!Point3D class methodsFon 'class initialization*! 

initialize 

RandomSouice «■ Random new 

"Point3D initialize"! ! 

!Point3D class methodsFon 'instance creation' ! 

x: xValue y: y Value z: zValue 

"Answer an instance of me with coordinates xValue, y Value, and zValue.' 

t self new setX: xValue setY: y Value setZ: zValue! 1 

PoinGD initialize! 



59 



EP 0 416 568 A2 



'From Smalltalk-80, Version 2.3 of 13 June 1988 on 21 August 1989 at 11:55:42 am'! 

Line subclass: #Line3D 

instance VariableNames: ff 
class VariableNames: M 
pooiDictionaries: M 
category: 'Graphics-Paths'! 



!Line3D methodsFor 'displaying*! 

dispLayOn: aDisplayMedium at: aPoint clippingBox: clipRect rule: anlnteger mask: aForm 
M The form associated with this Path will be displayed, according 
to one of the sixteen functions of two logical variables (rule), at 
each point on the Line. Also the source form will be first ANDed 
with aForm as a mask. Does not effect the state of the Path." 

collectionOfPoints size < 2 ifTrue: [self error 'a line must have two points']. 
aDisplayMedium 

drawLine: self form 

from: (self beginPoint) as2DPoinc + aPoint 
to: (self endPoint) as2DPoint + aPoint 
clippingBox: clipRect 
rule: anlnteger 
mask: aForm! ! 

!Line3D methodsFor: 'comparing'! 

<= aLine 

"compare the starting points of a 3D line" 

t self beginPoint <= aLine beginPoint! ! 
!Line3D methodsFor: 'converting'! 
as2DLine 

"convert 3D line to 2D line by romoving z coordinates" 

Itx ty aLinel 
aLine «- Line new. 
tx «- self beginPoint x. 
ty *- self beginPoint y. 
aLine beginPoint: tx@ty. 
tx *• self endPoint x. 
ty «• self endPoint y. 
aLine endPoint: tx@ty. 
taLinei 

asVector 

"convert 3D line to a vector" 

Ix y z aVectorl 

x (self endPoint x) - (self beginPoint x). 
y «- (self endPoint y) - (self beginPoint y). 
z «■ (self endPoint z) - (self beginPoint z). 
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aVecior «- Vector x;xy:yz:z w: 1.0. 
t aVector! ! 

!Line3D methodsFor 'transforming*! ! 

ILineSD methodsFor: 'testing'! 

sameAs: aLine3D 

"see if lines are the same" 

(self beginPoint = aLine3D beginPoint) 
ifTrue:[(seif endPoint = aLineSD endPoint) 

ifTrue:[ t true].], 
t false! ! 

it _ rtj 

Line3D class 

instance VariabieNames: "! 



!Line3D class methodsFor: 'examples'! 
sampleLine 

"A straight path with a square black form will be displayed connecting the 
two selected points." 

"Line3D sampleLine." 

t aLine aForm I 

aForm «- Form new extent: 5@5. "make a form one quarter of inch square" 

aForm black. "turn it black" 

aLine <- Line3D new. 

aLine form: aForm. "use the black form for display" 

aLine beginPoint: 100@100<S>25. 

aForm displayOn: Display at: (aLine beginPoint) as2DPoinL 
aLine endPoint: 600@600@100. 

aLine displayOn: Display. "display the line"! ! 

!Line3D class methodsFor: 'instance creation'! 
new 

"Answer a new instance of the receiver that is essential a single point 
at the upper left comer of the screen in 3D space." 

1 newSelf I 

newSelf *■ super new: 2. 
newSelf add: 0.0@0.0@0.0. 
newSelf add: 0.0@0.0@0.0. 
tnewSelf! ! 
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Point3D subclass: #Vector 

instanceVariableNames: 'w ' 

classVariableNames: " 

poolDictionaries: " 

category: 'Graphics-Circuits'! 
Vector comment: 

'This class are vectors for manipulating surfaces.*! 



[Vector methodsFon 'mathematical functions'! 

angle: a Vector 

"find the angle between two vectors." 

I theta c i 

c «• self cos: aVector. 

t theta (c arcCos) / 0.01745329251 



cos: aVector 

"find the cosine of the angle between two vectors." 



I c d ml m2 ) 

d <- self dotProduct: aVector. 
ml *- self magnitude. 
ml «- aVector magnitude. 
tc«-d/(ml*m2)! 



projection: aVector 

"Find the projection of the second vector onto the first vector." 



I aNew Vector dl d2 1 
dl «■ (self dotProduct: aVector). 
d2 «- (self dotProduct: self)* 
t aNew Vector «- self * (dl/d2)! 



sin: aVector 

"find the sine of the angle between two vectors." 



Icdml m2l 

d «- (self xProducn aVector) magnitude, 
ml «- self magnitude. 
m2 «- aVector magnitude. 
tc«-d/(ml*m2)! 



unit 

"find the unit vector of a vector." 



I aVector m 1 

m <- self magnitude. 

t aVector «- self /m! ! 



! Vector methodsFor: 'arithmetic'! 



* aScalar 

"Find the vector which is the product of a vector and a scalar." 
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I aVector I 

aVector «■ Vector x: (self x * aScalar) 

y: (self y * aScalar) 
z: (self 2 * aScalar) 
w: (self w* aScalar). 

t aVector! 
/ aScalar 

"Find the vector which is the quotient of a vector and a scalar." 
i aVector I 

aVector «- Vector x: (self x / aScalar) 

y: (self y/ aScalar) 
z: (self zl aScalar) 
w: (self w / aScalar). 

t aVector! 

dotProduct: aVector 

"Find the dot product of two vectors. The result is a scalar." 

label 

a* self x* aVector x. 
b «- self y * aVector y. 
c«-self z* aVector z. 
t a «- a +b -k! 

magnitude 

"Find the magnitude of a vector. The result is a positive scalar value." 
Iml 

m «■ self dotProduct: self, 
t m sqrt! 

xProduct aVector 

"Find the cross product of two vectors. The result is the normal vector 
of the two." 

lijkabcdef I 
a «• self y * aVector z, 
b«-selfz* aVector x. 
c«- self x* aVector y. 
a «- a - (self z * aVector y). 
b «■ b - (self x * aVector z). 
c «- c - (self y * aVector x). 
t Vector x: a y: b z: c w: 0.0! ! 

! Vector methodsFon 'accessing*! 

w 

"return w value for homogeneous representation" 
t w! 
w: aHoat 

"save w value for homogeneous representation" 
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t w «- aRoat! 

x 

"return x value for homogeneous representation" 
tx! 
x: aFloat 

"save x value for homogeneous representation" 
tx»- aFloat! 

y 

"return y value for homogeneous representation" 
ty! 
♦y: aFloat 

"save y value for homogeneous representation" 
t y *■ aFloat! 

z 

"return z value for homogeneous representation" 
t z! 
z: aFloat 

"save z value for homogeneous representation" 

tz*- aFloat! ! 

IVector methodsFor: printing 7 ! 

printOn: aStream 

"The receiver prints on aStream in terms of infix notation.' 

xprintOn: aStream. 
aStream nextPut: S@. 
y printOn: aStream. 
aStream nextPut: S@. 
z printOn: aStream. 
aStream nextPut: S@, 
w printOn: aStream! ! 

! Vector methodsFor 'converting'! 

asMatrix 

"Answer with a matrix version of self' 
I mat I 

mat «■ Matrix new: 1 by: 4. 
matatPoint: 1@1 put: self x. 
mat atPoint: 1@2 put: self y. 
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mat atPoint: 1@3 put: self r 
mat atPoint: 1<2>4 put: self * 
t mat! ! 

tt M| 

Vector class 

instance VariabieNames: ,f ! 



! Vector class methodsFor: 'Insxzzct creation'! 

x: x Value y: y Value z: zValue u Value 
"create a vector" 

I aVector I 
a Vector «- self new. 
aVector x: xValue. 
aVector y: y Value. 
aVector z: z Value. 
aVector w: wValue. 
r aVector! ! 
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'From SmaIltalk-80, Version 2.3 of 13 June 1988 on 1 September 1989 at 4:01:34 pirT! 



IMairix methodsFor 'converting'! 

asPoint3D 

"convert a matrix to a 3D point" 

I point f 

point *• Point3D new. 
point x: (self atPoint: 1@1). 
point y: (self atPoint: 1@2). 
point z: (self atPoint: 1@3). 
t point! 

asVector 

"convert a matrix to a vector" 

I vl 

v «- Vector x:(self atPoint: 1@1) 

y: (self atPoint: 1@2) 
z: (self atPoint: 1@3) 
w: (self atPoint: 1<S><1). 

t v! ! 
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